The 


ilosophical Magazine 


FIRST PUBLISHED IN 1798 
A Journal of Theoretical - 


Experimental and Applied Physics 


Vol. 4 June 1959 No. 42 
Eighth Series 


£1 5s. Od., plus postage 
Annual Subscription £13 10s. Od., payable in advance 


Printed and Published by 


TAYLOR & FRANCIS LTD 
RED LION COURT, FLEET STREET, LONDON, E.C.4 


THE PHILOSOPHICAL MAGAZINE 


Editor 
Professor N. F. Mort, M.A., D.Sc., F.R.S. 


Editorial Board 


Sir LAwRENCE Bracco, O.B.E., M.C., M.A., D.Sc., F.R.S. 
Sir Grorce Tomson, M.A., D.Sc., F.R.S. 
Professor A. M. Tynpatu, C.B.E., D.Sc., F.R.S. 


AuruHors wishing to submit papers for publication in the Journal should 
send manuscripts directly to the Publishers. 


Manuscripts should be typed in double spacing on one side of quarto 
(8x10 in.) paper, and authors are urged to aim at absolute clarity of 
meaning and an attractive presentation of their texts. 


References should be listed at the end in alphabetical order of authors and 
should be cited in the text in terms of author’s name and date. Diagrams 
should normally be in Indian ink on white card, with lettering in soft pencil, 
the captions being typed on a separate sheet. 


A leaflet giving detailed instructions to authors on the preparation of papers 
is available on request from the Publishers. 


Authors are entitled to receive 25 offprints of a paper in the Journal free 
of charge, and additional offprints can be obtained from the Publishers. 


The Philosophical Magazine and its companion journal, Advances in Physics, 
will accept papers for publication in experimental and theoretical physics. 
The Philosophical Magazine publishes contributions describing new results, 
letters to the editor and book reviews. Advances in Physics publishes articles 
surveying the present state of knowledge in any branch of the science in which 
recent progress has been made. The editors welcome contributions from 
overseas as well as from the United Kingdom, and papers may be published 
in English, French and German. ‘ 


Persea) 


The Growth of Grain-boundary Voids Under Stresst 
By D. Hutt and D. E. Romer 


Atomic Energy Research Establishment, Harwell 
[Received January 1, 1959] 


ABSTRACT 


The effect of combined hydrostatic pressure, P, and of uniaxial tension, 
o, on the rupture time of polycrystalline copper wire in the temperature 
range 400° to 500°c has been determined. All the specimens broke by 
intergranular fracture, due to the growth of voids along grain boundaries. 
Preliminary experiments indicated that the voids grow by the addition of 
vacancies under the action of the applied stress. A theory has been 
developed assuming that failure results from the growth of void nuclei on 
the grain boundary, which requires that the activation energy for failure is 
that of grain boundary diffusion and that the rupture time, t,, depends only 
on (o—FP). The experiments show that the activation energy is close to the 
expected value, and provided that o is constant t, agrees with theory. 
Changes in o, when (c— P) is constant, affect t, appreciably, and the results 
suggest that many of the void nuclei are stress induced. 


§ 1. INTRODUCTION 


In recent years attention has been drawn to the mechanism of failure of 
metals under certain creep conditions. Small voids have been observed 
at grain boundaries, particularly those transverse to the applied stress 
when specimens are tested at high temperatures ; rupture results from the 
growth and coalescence of these voids (Greenwood et al. 1954). This 
type of failure has been observed in Magnox cans used for containing 
uranium fuel elements, and the present work was undertaken to try to 
understand the mechanism of void nucleation and growth. 

There appear to be several methods whereby voids might be nucleated 
on grain boundaries. Grain boundary sliding can open up a hole at the 
junction of three boundaries or a ledge on a single boundary. Impurity 
particles which lack cohesion with the lattice may act as pre-existing voids. 
A dislocation pile-up breaking through to the boundary might also be a 
suitable nucleus. Of these four possibilities, all except the first must have 
vacancies supplied to them if they are to grow. Balluffi and Seigle (1957) 
have shown that sufficient vacancies can be produced at transverse grain 
boundaries because of the tension acting across the boundary, and suggest 
that this method of void growth is more likely to occur than that based 
on condensation of vacancies produced during strain (Machlin 1956), 
If the latter mechanism is operative voids would be produced by both 
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tensile and compressive strains. However, no vacancies can be generated 
at a grain boundary by applying a compressive stress across it. It was 
decided, therefore, to carry out a series of short time creep tests in 
which a uniaxial tensile stress and a hydrostatic pressure were applied 
simultaneously. 

In the first part of this paper we outline a mechanism of rupture based on 
the growth of void nuclei on a grain boundary by the addition of vacancies 
from the boundary. The effect of hydrostatic pressure and uniaxial 
tensile stress on the time for failure of such a system is calculated. The 
results of experiments designed to test this model are given and these are 
discussed. 


§ 2. THrory oF RupruRE BY VOID GROWTH 


The initial experimental results suggested that the mechanism proposed. 
by Balluffi and Seigle (1957) was probably the operative one for vacancy 
formation. This process must certainly be occurring, regardless of any 
other mode of void growth which may also operate, and it permits an 
analysis in terms of parameters which may be estimated with reasonable 
precision. 

In this section we shall derive an approximate expression for the rupture 
time, in order to show the significant features of this model. A fairly 
rigorous solution for a system having a simple geometry is given in the 
Appendix. 


2.1. Hstimate of Rupture Time 


The specimen is taken to be in the form of a wire and subjected to a 
uniaxial stress, o, along its length and a hydrostatic pressure P. Void 
nuclei of an unspecified nature are assumed to lie on the grain boundaries. 
Po will denote the radius of a typical nucleus and a their mean separation. 

At the temperature at which the tests are being made (400° to 500°c) 
grain-boundary diffusion is more efficient than lattice diffusion in trans- 
porting the surface atoms of the void to the grain boundary. This may 
be shown briefly by the following calculation : 

The ratio of the number of atoms transferred from the void to the 
boundary by lattice diffusion to that by grain-boundary diffusion is given 
approximately by 


D,.p 
1D eno! 
where D, and D, are respectively lattice and grain-boundary diffusion 
coefficients, p is the void radius and 62z the grain-boundary width. 
D,.6z~3 x 10-¥ cm sec for silver at 500°c (Hoffman and Turnbull 
1951), and a similar value is to be expected for copper. At the same 
temperature D, (for copper) ~ 2x 10-14 cm? sec. Even if p= Osc 
which must be an extreme upper limit to its possible size, the contribution 
from lattice diffusion is only 6% of the total, so that we are justified in 
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neglecting it. If there is considerable enhancement of lattice diffusion 
under stress, this argument may no longer be valid; the results, however, 
seem to indicate that lattice diffusion does not occur. 

In view of this it is presumed that the atoms are released from within a 
ring round each void, which lies in, and has the same thickness, 5z, as the 
grain boundary. 

The rate of growth of a void is determined by the gradient of chemical 
potential, V/, in the plane of the grain boundary, since the diffusion flux, /, 
is given by (Herring 1950) 
peas D, 
= rg 
where &=Boltzmann’s constant, 7’=absolute temperature, Q=atomic 
volume, and Vf is to be evaluated at the void surface. 

In the grain boundary, the potential is given by (Herring 1950) 


f= —-o,0 


where o, is the normal tension acting across the boundary. At a point 
well away from any void it may be assumed that this is roughly equal to 
the resultant of the externally applied stresses and has a maximum value 
(o—P) on a grain boundary normal to the direction of o. 

On the surface of a void f= — 2yQ/p, where y is the surface tension of the 
metal. 

Hence, as a rough approximation we may take, 


j 


giving 


; ID 2y 
Feta Bi Nola OPS 
: ina” =) 


This shows that only those nuclei for which 
2y 
Po> ied P 
will grow. For voids smaller than this critical size the applied stress is not 
sufficient to overcome the surface tension and sintering occurs. Assuming 
that there are some nuclei which can grow, the larger ones will determine 
the rupture time and it then is reasonable to neglect 2y/p with respect to 
(o—P). 
The area of the void surface from which diffusion takes place is 2zpéz, 
so that the rate of increase in void volume is given by 
dV (D,8z)(o—P) + 27pQ 
dt kTa 
If the bubble retains spherical shape by surface diffusion, we have 
dp | (D,2z)(o — P)Q 
dt 2kT ap ; 


2X2 
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To obtain the rupture time, ¢,, this expression must be integrated between 
p=py and p~ 3a. 
Assuming that a> py, this gives 
kT a8 
r™ 4(D,o2)(0 — P)Q." 
It is shown in the Appendix that for an array of equally sized nuclei 
lying on a square lattice the rupture time is given almost exactly by 
ae kT a 5 Leeann 4,2 ia | ea ~ Al 
 27Q(D,62)(o—P) (8 (¢-—P) (e—P}* pola —P) — 2y 


t 


(11) 


The approximate result differs by a factor of 47 from the first term in the 
bracket in eqn. (11). The effect of nucleus size, which was neglected 
above, appears in the third term and is usually negligible except for values 
of py only slightly greater than 2y/(o—P). 

It would appear that the square array is about the best simple periodic 
approximation that can be made to a random arrangement of nuclei, 
but in order to estimate the effect of the assumed spatial distribution on 
the result it was decided to carry out an analysis similar to that given in 
the Appendix, but for a single void feeding atoms into a concentric ring 
of grain-boundary having the same initial area, 1.e. (a2—7p )”). A close- 
packed array of such systems is an alternative approximation that may 
be made to a real grain boundary. Using a particular set of parameters 
the time for failure on this model was about 70% of the breaking time 
using the square lattice. The longer time taken by the square lattice 
is due to the fact that the region of grain boundary associated with each 
void is of a more irregular shape than that of the other model (in which 
circular symmetry is used), and is therefore the more realistic. 

It therefore seems reasonable to use eqn. (11) for the interpretation of 
experimental results, with the following proviso. A randomly dispersed 
array of mean separation a should bring about fracture in a shorter time 
than the uniform set considered in obtaining eqn. (11), as some will join 
up quickly and increase the local stress concentration. A value of a 
determined from an experimental result will therefore be somewhat less 
than the mean distance between nuclei. 


2.2. Theoretical Predictions 
The significant features of eqn. (11) are: (i) The activation energy of 
the process is that for grain-boundary diffusion. (ii) The quantities o 
and P only occur in the form (o—P) so that for identical specimens at the 
same temperature the breaking-time is a function only of (s—P). If 
(o—P)> 2y/pp the three terms on the right-hand side of eqn. (11) are in the 
approximate ratio 


1 2 
ot ee | A 
8 a(o—P) a(a—P) 

Taking y=10?dynecm-!, a~2x10-4cm, the first term dominates 
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if (o—P)25x 107 dyne em-*=750 p.s.i. In the experiments to be 
described this is true, although significant contributions come from the 
other two terms. Thus the shape of the t, versus (o—P) curve is 
approximately hyperbolic down to (o—P)= 2y/p9, at which point the third 
term becomes infinite. Physically this is the point at which the applied 
stress is balanced by the surface tension of the void nucleus, so that the 
voids never grow. Figure 1 shows the expected variation of t, with (o— P). 


Fig. 1 
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28 ipo (6-P) 


Theoretical Variation of t, with (o—P). 


§ 3. EXPERIMENTAL DETAILS 


The tests were made on a relatively impure copper (99-8°%) which 
contained small copper oxide inclusions in short stringers along the direction 
of drawing. The specimens were in the form of wire 0-020 in. diameter, 
which was etched in dilute nitric acid to reduce the diameter to between 
0-013 and 0-010 in. along a 1 in. gauge length. Copper beads were welded 
onto the ends of the wire specimens to hold them in the apparatus. The 
specimens were annealed at 600°c for one hour before testing. Combined 
hydrostatic and uniaxial tensile stresses were applied using the apparatus 
illustrated in fig. 2. The tensile stresses were obtained using a spring D 
attached to one end of the specimen C and held in a silica tube B, approxi- 
mately 12 in. long. This assembly was inserted inside a pressure bomb A 
(Fennell et al. 1958) with an internal diameter 5/16 in., capable of holding 
pressures of 6000 p.s.1. for a few hours. All the tests were made in argon. 
The time to fracture was measured by means of an electric clock circuit 


678 D. Hull and D. E. Rimmer on the 


through the specimen. The pressure bomb was inserted in a furnace tube 
and temperatures between 400° and 500°c maintained to +3°c. The 
creep strain was not measured during the test, but measurements on 
fractured specimens showed that the creep strain was less than 3 to 5%. 
The specimens were examined metallographically after fracture by 
sectioning along their length followed by polishing on fine alumina. 


Diagram of apparatus used for applying simultaneous hydrostatic pressure 
and uniaxial tension. 


§ 4. EXPERIMENTAL OBSERVATIONS 


4.1. Effect of Non-hydrostatic Compressive Stresses 


An experiment was carried out to determine whether voids were formed 
in specimens tested with c=P. In this case there is no resolved stress 
along the tensile axis, but there are compressive stresses normal to this 
axis. The specimens were tested at 480°C with a tensile stress sufficient 
to break the specimens in less than 20 min (without the hydrostatic 
pressure). The specimens did not fracture after 8 hours and subsequent 
micro-examination revealed that no voids had formed. 


4.2. Metallography 

It was confirmed that all the specimens tested had ruptured by the 
formation of voids on grain boundaries. Figure 3 (Pl. 77) shows the 
general appearance of a specimen, the majority of voids having formed on 
boundaries normal to the applied stress. A few stringers of oxide particles 
are also visible. Figure 4 (Pl. 77) is a high magnification photograph 
which shows clearly the voids along grain boundaries. Voids do not form 
on twin boundaries, even when the boundaries are transverse to the 
applied stress. The photographs show that the voids in fractured specimens 
are close together and often coalesce. The size and number of voids 
varied with the applied stress. This is illustrated in the photographs 
(figs. 5 and 6 (Pl. 78)) of two specimens tested at the same temperature, 
the first stressed at 4730 p.s.i. and the second at 3000 p.s.i. At the high 
stress a large number of small voids have formed compared with a smaller 
number of larger voids at the low stress. 
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4.3. Activation Energy of the Process Leading to Rupture 


A series of specimens was tested at a constant uniaxial tension and the 
time to fracture measured for a range of testing temperatures between 
370° and 490°c, The results, which are plotted in fig. 7, give an activation 


Fig. 7 
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Variation of log ¢, with 1/7 for c=4730 p.s.i., P=0. 


energy of 25 000+ 3000 cals/mole. This figure compares well with the 
activation energy for grain-boundary diffusion determined for silver 
(20 200 cals/mole) (Hoffman and Turnbull 1951). 
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Fig. 9 
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Variation of ¢, with o for constant values of (o— P). 
(Taken from smoothed experimental results at 410°c ) 
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4.4. Effect of Stress (oc) and Pressure (P) at Constant Temperature 


The variation of the time to rupture (¢,) with hydrostatic pressure for 
a number of values of o was obtained at 410°c and the results are plotted 
in fig. 8. The variation of rupture time with stress for a number of values 
of (c—P) has been obtained from these results and is plotted in fig. 9. 
Equation (11) predicts that the rupture time should be dependent only 
on (o—P), i.e. a variation of ¢ should not affect the rupture time providing 
(o—P) remains constant. As shown in fig. 9, the rupture time alters 
appreciably with change in o and this is discussed below in § 5. 


§ 5. Discussion 


In $4.1 an experiment designed to test the theory of void growth by 
vacancy condensation under strain was described. It was shown that 
although a certain stress applied in tension would easily break the specimen, 
when applied in compression it could not doso. Further, in the latter case 
no visible voids were produced. This suggests that the vacancies for 
void growth are produced by the tension across the grain boundaries 
rather than deformation inside the grain. It seemed reasonable therefore 
to try and correlate the experimental results with a theory based on the 
grain-boundary diffusion model from which eqn. (11) was derived. 

The activation energy of the processes leading to rupture in the tem- 
perature range used corresponds with that for grain-boundary diffusion 
rather than lattice diffusion. This supports one of the basic assumptions 
of the theory that the voids grow by the diffusion of atoms from the void 
surface along the grain boundary. 

Examination of the results in fig. 8 compared with eqn. (11) shows that 
the theory cannot explain all the experimental results, since if (o— P) were 
the only stress factor affecting ¢, then for a constant value of (o — P), e.g. 
3000 p.s.i., the breaking time should be constant, whereas at o = 4730 p.s.1. 
t.=2-2x 104 sec, and at o=5050 p.s.i., t,=1-25 x 10* sec. However, we 
shall consider first how eqn. (11) agrees with the effect of (c—P) on #, 
providing o is constant, and then examine the variation of ¢, with o. 
Figure 8 shows some observed values of ¢, as a function of P for a number 
of values of c. Equation (11), with which it is to be compared, contains 
two unknowns, @ and py. The appropriate value to use for py is difficult 
to determine. It is however possible to find a lower limit for it because it is 
known that tf, 00 as (o— P)->2y/py (see § 2.2). Figure 8 shows that ¢, is 
still finite when (c— P)=2000 p.s.i. It is unlikely that this, the smallest 
experimental value, is by coincidence very close to the critical one. If we 
take 1500 p.s.i.= 108 dyne cm as the critical value 


2x 103 


a = —5 
5 “He J On Gm 


which isreasonable. Using this figure, any one of the experimental results 
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giving t, for specified values of o and P will enable a value for a to be calcula- 
ted by eqn. (11). This may then be used to predict values of ¢, for compari- 
son with the remaining experimental values. If results for which P=0 
are used to determine a, the term containing py is small and the assumptions 
made in its estimation have little effect on the calculated value of a. 

It was found that for each value of c, a different value of a is needed 
to fit the experimental point at P=0. However, if for any given a, the 
appropriate value of a is used to derive the ¢,(P) relation, this agrees well 
with the experimental points. The full lines in fig. 8 are the theoretical 
curves, using the following values of a: 


Table Giving the Distance Between Voids, a, for 
Values of Applied Tension, o 


As explained in § 2.1 the values obtained for a are likely to be rather less 
than the mean distance between nuclei. This suggests a void separation 
of about 3 y, a figure which is in general agreement with the metallographic 
observations (see figs. 3 and 4). 

The effect of variation of o on t, for constant values of (o—P) is shown 
in fig. 9. At higher values of o the rupture time is less. This is in agree- 
ment with the results given in the table, where it is shown that the distance 
between voids a, decreases with increasing o and therefore if the growing 
force (c—P) is constant the rupture time will decrease with increasing o. 
The observation (§ 4.2) that there are more voids present after a tensile 
stress of 4730 p.s.i. has been applied than after 3000 p.s.i. supports this 
conclusion. 

Thus it would appear that the mechanism of void growth is as postulated 
in the theory, but that in order to explain all the experimental facts it is 
necessary to allow a density of void nuclei which increases with applied 
tension. ‘This virtually implies that the stress itself is generating nuclei. 
These could result from grain boundary sliding as suggested by Gifkins 
(1956) and Chen and Machlin (1956). It is also possible that some of the 
impurity particles are more firmly bound to the lattice than others, and 
only become free when the specimen is strained. 

It does not seem possible to make any useful deductions about the 
number of nuclei produced by stress. Estimates made from the a values 
given in the table would require the assumption that all nuclei are produced 
at the beginning of the experiment. This also prevents a useful comparison 
with the observed variation of void separation with stress. 
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In conclusion, the experiments show that the grain-boundary voids 
which are frequently observed to be the cause of creep failure, grow by the 
diffusion of surface atoms to the surrounding grain boundary under the 
action of the applied stress. It appears that many of the void nuclei are 
induced by the stress, but their exact nature is uncertain. 
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APP iN DIX 
The Rupture Time for a Square Array of Equally Sized Nuclei 


Figure 10 shows a square array of spherical voids with repeat distance a 
and having grown to a radius p, lying on a grain-boundary normal to the 
direction of c. The centre of a void is taken as the origin of rectangular 
and polar axes. 

In order to calculate Vf at the void surface it is necessary to define the 
problem more exactly than in § 2, and we shall make the following additional 
assumptions : 

(i) The ring of void surface lying in the grain boundary, which acts 
as a source of diffusing atoms has a radial width w, about equal to that of 
one atomic layer (see fig. 10). The atoms are released uniformly from 
within this volume. 

(ii) The atoms are deposited uniformly over the grain boundary, so 
as to minimize lattice strains. . 

(iii) The void retains its spherical shape throughout its growth. It 
must be pointed out that whiist there is mechanism available (i.e. surface 
diffusion) for doing this, a sphere is not the equilibrium shape for a void 
under the conditions that have been postulated. Owing to the presence 
of the other voids the natural shape of its grain-boundary cross section 
is a deformed circle with bulges towards the corners of the unit square. 
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These become more predominant as the void grows. In order to retain 
a circular cross section it is necessary to impose small restraining forces 
onthe void. In view of the other approximations necessary in this problem 
it has been thought permissible to assume that the bubbles remain spherical, 
yet ignore the forces holding them to this shape. Although the effect 
of this on the final result is unimportant, it will be seen later that it produces 
some slight inconsistencies in the mathematics which have been accepted 
in view of the great simplification in the working that has resulted. 


Fig. 10 


x 
Cross section through voids lying in the plane of the grain boundary. 


Let us consider a unit cell centred on the origin and defined by 
|x| < $a, 
ly| <a. 
V?f is a measure of the strength of a given point as a source of or sink for 
atoms. In view of assumption (i) and (ii) above it is determined (apart 
from a scaling factor) at all points inside the unit cell since we have, 
V3f=0 when 0<r<p, ] 
Vij= A when <7 <p--40.. +) se een) 
V4f=B when r>p+w, | 
where A and B are related by the equation, 
B(a? — mp?) + Aw: 2ap=0 ae = Se ee) 
as the total number of atoms is constant. 
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But V?f has periodicity a in a, y, and may be written in the form 


Vif= Sy S Cyn C08" (mee +ny) eee Sat) 


m=—on=—0o 


where 


+ia ¢+ta On 
| Vf cos — = ne + ny) dx dy 


—ta 


B ftie on a 
=a] “ ere, ae cos “a 


(m cos 8+ n sin 0) dé me 
A p+w (27 by 
+= | cos (m cos +n sin 9) dé rdr, 
a Jo 0 a 
from eqns. (1). 


Assuming w <p, 
even & a pB 2mp(m? + n?)12 
Cmn = (B or 0) alm? nie Jy ia | a 
2 > 2 2\1/2 
+ a ie — ) io Sena} 
where the first term is B if m=n=0 and zero otherwise. 
Hence 


5) 2 
Omg= B— OR AE 4 Smee =B(1- LE yes 
a~ a a 


a a 
from eqn. (2). Otherwise, 
Con ete p di Qarp(m? + n2)H2 fae al J Qmrp(m? + n2)U2 
B a(m?+n?)H2 > + a Arn A eueeare eames) 


(5) 


The most te ee of eqn. us ) which has the required symmetry is 


2Qar 
K ne ee 
pad Dati 08S (mat ny) + (6) 
where K is a constant, sso from the condition, 
f=—¢,2. Re he ets eee ) 


o,, Will vary over the unit cell owing to the stress concentration set up round 
the void. However, the total force over the unit cell due to the applied 
stress and pressure must be (o— P)a?. 

Hence integrating eqn. (7) over a unit cell and allowing for surface 


tension, 
+ia p+ta p (2a 
Oe (e=P)a l= | ae ay" i Vari ee 1) 
—taJ —tha 
On physical grounds we may expect f to have the value —2yQ/p at all 
points inside the void. Ifit were calculated from eqn. (6) 6) this result would 
not be obtained because the restoring forces round the void have been 
neglected. —(2yQ/p) will in fact be the mean value inside the region and 
for the purposes of carrying out the second integration in eqn. (8) we shall 
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take it as constant. The origin is one of the positions which has the mean 
value, a fact which will be used later. Hence 


Q[27py — (o — P)a?] = Ka? + 20 . 7p” 
p 


giving 
K=—-—QO(c—P). 


Hence 2 
= - Slots ata os — . (ma+ny)—-Q(o—P). . . (9) 
Assuming a radial Hie ‘of atoms at the void surface which is reasonable 
if p <a the diffusion flux is given eg 


Lp of 
ie ah OF) pee 


If we let w+0 and evaluate (f/dr) at r=p it must be remembered that 
(Af/0r) is now discontinuous at this point (see fig. 11) and the value required 
is the maximum (numerically) which is equal to twice that obtained by 
computing the Fourier sum at r=p. 


Fig. 11 
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Approximate radial variation of gradient of chemical potential. 
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The number of atoms, per second, coming out of an elementary area 
pdész of void is, 


_ 2D, pd0dsz of 
are Or) 
Therefore, the rate at which atoms leave the void is given by, 


aN _ _ 2p(D,82) 52) or of do 
dt TT QS 0 or r=p 


ar’ 2(D,82)pa Cran oy. (a xr oe) (10) 
RPQ 2% (mB nye a | 
The coefficients C),,, are not known, but ©,,,,/B is given in eqn. (5). 
However, as stated earlier the value of f is known at the centre of the void. 
Substituting 


f=- ie when x=y=0 
Pp 
into eqn. (9) gives 
aN _ pt = (0-P- x) aN 
P 


dt kT a 
where 


B 2a7 

ee Crnl p 2 2) =a 
a NT ms & aed Z| 225 go 
The rate of change of void radius is, 

dp _  Q (dN _ 2nQ(D,62)(o—P —2y/p)S 


dt 4p? dt kT ap 
The time to fracture is obtained by integrating this expression between 
p=p, and p=a. S isa function of p/a only, and has been calculated for 
a set of values in the expected range of integration. These suggest that 
the approximation S=1 may be used as it only deviates significantly from 
this value for a short time just before the voids coalesce. 


This gives, 
a a 62) “a= A 2 pep 
,o—P—2y/p 


where t, = time to ce a 
On integrating and assuming py <a@, we obtain 
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The Lattice Thermal Conductivity of Dilute Alloys of 
Silver and Goldt 


By G. K. Wurrrt, 8. B. Woops and M. T. Eirrorp§ 
[Received January 6, 1959] 


ABSTRACT 


The thermal conductivity of certain alloys was measured on samples 
which had been carefully annealed to reduce the concentration of disloca- 
tions. After subtraction of an electronic component estimated from the 
residual electrical resistivity, the remainder was proportional to 7? below 
10°x, and agreed with Makinson’s calculation of the lattice conductivity, 
assuming that lattice waves of all three polarizations interact equally with 
conduction, electrons. 


THE lattice thermal conductivity, A,, of metals at low temperatures appears 
to be limited by electron scattering of the lattice waves or phonons and 
measurements of A, may therefore yield information about the electron— 
phonon interaction (Klemens 1954a, 1956). The lattice conductivity of 
pure metals cannot be measured directly because of the dominance of X,, the 
electronic component of thermal conductivity, but it can be deduced under 
certain assumptions from measurements on dilute alloys in which the two 
conductivity components may be comparable. Studies, below about 
100°xK, of the electrical and thermal conductivities of silver alloys (Kemp et 
al. 1954, 1956), of copper alloys (White and Woods 1954, Kemp, Klemens 
and Tainsh 1957, Kemp, Klemens, Tainsh and White, 1957) and of a gold 
alloy (Birch et al. 1958) have yielded values for the lattice thermal con- 
ductivity of silver, copper and gold. These values lead to the conclusion 
that the interaction of the conduction electrons with transverse lattice 
waves is almost as strong, if not as strong, as with longitudinal waves. 
Unfortunately the interactions of lattice waves with dislocations and 
with conduction electrons produce components of thermal resistance which 
vary in the same way (oc7’*) with temperature. Thus, if there is an 
appreciable number of dislocations in the specimens, it is only possible to 
set an upper limit to the thermal resistance produced by the electrons. 
Although annealing at high temperature may substantially reduce the 
dislocation density, it is possible that many dislocations become immobile 
in alloys having large concentrations of solute atoms. However, in suffici- 
ently dilute alloys suitable annealing should remove all but a negligible 
number of dislocations. Some years ago measurements were made on a 
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dilute solution of iron in copper (White and Woods 1954) and this work has 
now been extended to dilute alloys of iron in gold and tin in silver. 

The apparatus and techniques described by White and Woods (1955) were 
used to measure the electrical and thermal conductivities of the specimens 
between 2 and 150°x. The alloys were vacuum cast from the constituent 
metals and homogenized in this laboratory and were then swagedt} to form 
rods which were finally sealed in evacuated quartz tubes and annealed. 
Leads were connected to the specimens with a bismuth-cadmium solder 
which does not become superconducting above 1°K. 


Table 1. Some Physical Data for the Specimens 


Solute Length | Diameter 


Specimen Anneal 


‘ P295 Po . WoT " 
(at. %) (em) (wQ em) | (uO em) | W-t em deg? 


Au | 0-13 % Fe 3°48 2: 15 hrs 700°c 3-44 “18, 
Ag l 0-14% Sn 3-10 2-4/ 15 hrs 700°c 2-32 “724 
Ag 2 0-3% Sn “f : 150 hrs 760°c 2-69 c 


Some physical data for the specimens are given in table 1. Te check the 
compositions, part of each alloy was chemically analysed after being 
annealed. Also the electrical resistance values measured at 295°K and in 
the liquid helium range were compared with the measurements by Linde 
(1939) on alloys to estimate the solute concentration (Gerritsen 1956). 

The total thermal conductivity of the alloys is given by the sum of the 
electronic and lattice components, 


A=A,+A, FO ie ee at ke Pee ES) 
where the electronic thermal resistivity is 
LA = W,= Wot W,. Sect a0 Seo, see Ce) 


The ideal electronic thermal resistivity, W,, which arises from electron 
scattering by lattice vibrations, is assumed to be the same as that of the 
pure parent metal. The residual electronic thermal resistivity, Wo, due 
to impurity scattering, should be related at temperature 7’ to the residual 
electrical resistance, po, by the relation 


AW ATL Meee, er ae(3) 


where the Lorenz number, L,= 2-45 x 10-8 watt ohm/deg”. At tempera- 
tures below about 20°K, W,<W y= 1/Ao, and thus A=A,)+A,. Since we 
expect \,=constant x 7? (see Makinson 1938, Klemens 1954 a) at low 
temperatures where the lattice waves are scattered mainly by electrons, A/T 
plotted versus 7’ should yield a straight line with intercept AG anes 0, 
The values of W,7' given in table 1 are obtained in this way and agree within 


+ We are indebted to Mr. J. A. Perry of the Mines Branch, Department of 
Mines and Technical Surveys, Government of Canada for arranging to have the 
alloys swaged in his department. 
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about 1°% with the values obtained from eqn. (3). Values for W, were 
taken from measurements on pure gold (White 1953 a, Rosenberg 1955) and 
on pure silver (White 1953 b, Rosenberg 1955). Using eqn. (2) we were then 
able to deduce A, and hence, with eqn. (1), A, The results are shown in 
fig. 1 for gold and fig. 2 for silver. 

It can be seen that below about 10°K the lattice conductivity does indeed 
vary as 7” and the constant A ,7'* derived from the results in this region may 
be compared with theoretical values. Klemens (1954a) has derived two 


Fig.l 
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values theoretically—one on Makinson’s (1938) assumption that the elec- 
trons interact equally with lattice waves of all three polarizations, the other 
adopting the Bloch model in which interactions occur with the longitudinal 
waves only. 

Makinson’s coupling scheme results in a lattice conductivity about twenty 
times higher than that given by Bloch coupling and, although an electron 
interaction with transversely polarized waves intermediate between these 
two extremes would lead to intermediate conductivity values, it is apparent 
from table 2 that the experiments support Makinson’s hypothesis very 
closely. This result generally confirms earlier work on less dilute alloys, the 
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present values of \,7'~* being only a few per cent greater than the maximum 
values obtained earlier of about 1-9 x 10-8 for Ag (+ 2% Cd) by Kemp et al. 
(1956) and a recent value of 2-6 x 10-3 for Au (+0-7°% Pt, annealed 750°c 
and 1050°c) by Birch, Kemp, Klemens and Tainsh (to be published). This 
supports the thesis that provided the solute content is less than of the order 
of 1% and the alloy is annealed under normal conditions, then insufficient 
dislocations are locked into the alloy to noticeably affect the measured 
lattice conductivity. Hence a value of \,7'—? (or W gl”) representative of 
phonon-electron interaction is obtained in these dilute alloys. 


Fig. 2 
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Near liquid air temperatures scattering of the lattice waves either by 
mutual interactions due to anharmonicity or by point imperfections should 
lead toA,c 7. The observed lattice conductivity exhibits this tempera- 
ture dependence but is smaller than theoretical estimates by a factor of four 
or five (see White and Woods (1958) for a discussion of anharmonicity 
theories; Klemens (1956) for point scattering theory). Although uncer- 
tainties in the theoretical estimates make it difficult to draw conclusions 


2Y 2 
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Table 2. Experimental and Theoretical Lattice Conductivity at Low 


Temperatures (< 10°K) 
62 Re ie a ee 


dg’? (watt em-1 deg-*) 


Metal 

Experimental Makinson Bloch 
Gold 2:8 x 10-4 2-9 x 10-3 Lb cl 0—4 
Silver 2-0 x 10-3 2-3 x 10-% 1-2 x 10-4 
Copper 1-8 x 10-3} 1-4 x 10-3 Ven 10 


+ See White and Woods 1954. 


from this comparison it may be worth noting that in the inert gas solids, 
which also have a cubic lattice structure, somewhat better agreement 
is obtained (White and Woods, 1948) with the anharmonic interaction 
theories. 
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ABSTRACT 
A completely theoretical calculation was made of the dependence of the 
yield point on the degree of long-range order in an AB b.c.c. super-lattice. 
The resulting equation is 
t= NkTS,2/4b 


where 7’ is temperature, and Sy is degree of long-range order. For B-brass 
@ maximum in 7 at 515°K of 0-54x 10° dynes/em? is predicted; the experi- 
mentally observed values are 493°K and 0:36 x 10° dynes/em?. 

A sharp increase in 7 was observed near 7',; it is associated with short- 
range order strengthening as predicted by Sumino. 

The width of an anti-phase domain boundary was also calculated. For 
T>0-5T; the width is given by the equation 


width = 4n[T./(To—T)]!2+4. 


§ 1. INTRODUCTION 


THE idea of the super-dislocation in a super-lattice was first presented 
by Koehler and Seitz in 1947. 

The effect of long-range order on the yield point was first discussed 
quantitatively by Cottrell (1954), who considered the size effect of an 
anti-phase domain. The Cottrell equation is 

mcs (1 bl 1) Ne ae ee er oP, ee GL) 
where 7 is yield stress, y is surface energy of a domain boundary, / is size 
of the domain, a is width of the boundary, and a~6 is a shape factor. 
The equation is based on the idea that the external work 7z/ equals the 
change in energy on the slip plane from yaal to y. Equation (1) predicts 
a maximum in the yield point at a domain size of about 20 4°. Implicit in 
eqn. (1) is the assumption that yielding is associated with a shear 
displacement of 1/2 along the slip plane. In general, it would be expected 
for large | that the strain would be so large after a shear of //2 that the 
resulting stress would be more closely associated with a flow stress rather 
than a yield point. Consequently, eqn. (1) should be restricted to small 
values of 1. 

A recent theory by Sumino (1958) considered two effects: (1) the 
disordering of short-range order as proposed by Fisher (1954) and (2) the 
interaction of the stress field of the dislocation with the change in nearest 
neighbour distance caused by ordering. The disordering of short-range 
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order is the larger of the two effects and its corresponding yield rises very 
rapidly as the temperature approaches 7’, from above. The other effect 
is such that screw dislocations are much easier to move than edges, and 
its corresponding yield stress is a maximum at 7’,. 


Fig. 1 
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Paired dislocations joined by an anti-phase boundary. 


Ardley and Cottrell (1953) proposed that in the super-lattice the stable 
dislocation consists of a pair of equal dislocations joined by an anti-phase 
boundary (fig. 1). The equilibrium configuration of such a super- 
dislocation was calculated by Brown and Herman (1956) for an AB b.c.e. 
super-lattice. Neglecting the entropy, the equilibrium spacing for two 
edge dislocations is 

r= 4Gb*/34/2n(1 — v)kTS,?. aR, Ree | 
When So, the degree of long-range order, is equal to unity, r becomes 
32 Burgers vectors for 6-brass, and the energy of the anti-phase domain 
boundary 82ergs/cm?. The stress to move a super-dislocation has been 
discussed many times but a concrete theory has not been proposed. The 
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general attitude indicated that the stress to move a super-dislocation 
should not be appreciably greater than for an ordinary one because the 
trailing dislocation should re-order the disorder produced by the leading 
one. In this paper the nature of the anti-phase domain boundary in the 
super-dislocation is investigated further; and it turns out that the stress 
required to move it may be very large. The new mechanism of yielding 
predicts a maximum in the yield point as a function of the degree of 
long-range order. It is interesting that all quantitative theories (Cottrell 
1954, Sumino 1958) involving the effect of long-range order on yield point 
predict a maximum. Experimental data are presented which support 


the new theory and which also indicate that the effect predicted by 
Sumino is observed near 7’. 


§ 2. THEORY 

It is a geometrical necessity that an anti-phase domain boundary 
starts or terminates at a dislocation. In the AB b.c.c. type of super-lattice 
every anti-phase boundary which starts at any $[111] dislocation can be 
terminated by any other }[111] dislocation. The same generalization 
does not hold for a f.c.c. super-lattice because it has four kinds of lattice 
site whereas the b.c.c. lattice has only two. The present theory will be 
confined to the AB b.c.c. lattice and this will be directly comparable with 


(a) (0) 
(a) Movement of paired dislocation not on same slip plane. AB is initial 
position; A’B’ is position after slip. 
(b) Movement of paired on same slip plane. AB is initial and A’B’ is position 
after slip. 


the Sumino theory. In general, dislocations not on the same plane can be 
joined by a curve anti-phase boundary; this case will not be considered 
because the stress to move these dislocations is much greater than for the 
case where the dislocations are in the same slip plane. In fig. 2(a) the 
super-dislocation consists of two edges on different slip planes, and when 
it slips, the anti-phase boundaries AA’ and BB" are generated. In 
accordance with Fisher’s equation, the required stress is 


ado ee Rae Le) 
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where H,,, is the surface energy of AA‘ or BB’. In fig. 2 (b) the super- 
dislocation lies on the same slip plane and the stress required for slip is 


T= (Epp — ‘AA’)/ 20 . 5 ° . . . 5 (4) 


where Epp and E,,, are the surface energies produced by the slip of the 
leading and trailing dislocation respectively. Consequently if, as it has 
been assumed in the past, #,),= 2,4, then no additional stress is required 
to move a super-dislocation. The main purpose of this paper is to prove 
that the stress to move a super-dislocation is in general not equal to zero. 
The proof requires that the anti-phase boundary generated by slip in the 
absence of diffusion is significantly different from an anti-phase boundary 
formed under thermal equilibrium (fig. 3). The boundary formed by slip 


Fig. 3 


Width of anti-phase boundary showing order versus distance from slip interface. 
Full line is anti-phase boundary generated by slip of a single dislocation. 
Dotted line is an anti-phase boundary formed under thermal equilibrium. 


is a step-function, but the one formed under thermal equilibrium is 
smoothly varying above 0°K as demanded by considerations of the 
entropy. Consequently, there is a difference in order on the slip plane of 
the two types of anti-phase boundaries. Since the degree of order 
determines the surface energies in eqn. (4), the stress to move a super- 
dislocation is in general not equal to zero except when the equilibrium 
order is unity or zero, Therefore, a maximum in the stress is expected 
at some intermediate degree of order. In order to determine the yield 
point we need to know the width of the anti-phase domain boundary 
formed under thermal equilibrium. 
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§ 3. CALCULATION oF WipTH oF ANTI-PHASE BOUNDARY 


A one-dimensional calculation will be made in which it is assumed that 
the order is constant along planes parallel to the slip plane. Obviously, 
the assumption breaks down at the core of the dislocation, and the 
assumption is more valid the more extended the super-dislocation. 
Equation (2) as applied to CuZn indicated that the minimum spacing at 
order S)= 1 was 320 ; therefore the one-dimensional calculation is expected 
to give results which have physical significance. Using the Brage—Williams 
parameter for long-range order, the fraction of A atoms on « sites of the 
nth plane from the slip interface is given by 


Considering only nearest neighbour interactions the energy per unit area 
of anti-phase boundary is given by 


EH=2N (2 SS 1 +3)@Vant Vaat Vpp) 


n=1 


sal > S,2+ SySye1—J8i?)(2Van—Vaa—Vaw)} (0) 


n=1 
where 2 is the density of atoms on the (110) plane, Vy,, Vgpg, and Vaz 
are the energy of the bonds between A—A, B—B, and A-B atoms; and 
S,, =p, the degree of long-range order in the perfect crystal under 
thermal equilibrium. In deriving eqn. (5) the widths of anti-phase 
domain boundary on both sides of the slip interface are included 
(S,,= —S_,). The entropy per unit area of anti-phase boundary is 


S=2k > log{NY3(1—S,)NI4(1+S8,)N. . . . (6) 


n=1 
The free energy is given by 
F=E-TS. eee ee os eee) 
Substituting eqns. (5) and (6) in (7) and setting dF'/dS,, =0, the following 
values of S,, at thermal equilibrium result : 
S74 Sos 27) 7 log (1 8,)/(1—8)) 


for n>1 Sy 2S ee ol log (la SL Sy) 2 ae. a(S) 
where kT .= —2(2Vap—Vaa— Vopp)- 
Equations (8) reduce to the well-known equation (Kittel 1953) 

So= 2/27 log (VF S5)/(1— So) ie Sur! ree 


where S, is the equilibrium degree of long-range order in the perfect 
crystal. The exact solution of eqns. (8) over the entire temperature range, 
0°x-7',, can only be obtained by numerical methods. However, since 
the higher temperatures are of the greatest interest, a high temperature 
approximation will be used which leads to a simple solution and which 
was found to agree closely with more exact calculations down to about 


0-57". 
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At high temperatures where S, <1 eqns. (8) transform to the differential 
equation 
48--atd?SidX*=47TS/T, . - « «© =~, (10) 
where S,, transforms to S(X), a is the distance between (110) planes, and 
X is the coordinate perpendicular to the slip interface. From fig. 3 the 
boundary conditions are S=0 at X=0 and ds/dX=0 at S=S,). A 
sinusoidal solution is indicated, 
S=S,sinpX; p?=4(T,-—T)/eT,. ie, eee 
The width of the anti-phase domain boundary is determined by the 
condition 
p(n ~—1/2)a=1/27 
and 2n.,, the number of atomic planes in the boundary, is 
20 5S 1f2a[T Eat a ee) 
n.. behaves properly because n~2 when 7’=0 and n,,>c when T+T’,. 
At 0-975 7, where S,=0-52n,,~10. Thus, the width of the anti-phase 
domain boundary is generally very narrow except very close to 7',. As 
far as the author knows the above calculation represents the first one for 
the width of an anti-phase boundary. A more accurate solution of eqns. (8) 
was made by approximating the region of S<1 with a sine function and 
the region of (S,—S)<1 with an exponential. The two solutions were 
matched at a slope of tanz/4. The increased accuracy over the purely 
sinusoidal approximation was found to be small and not worth the loss 
of simplicity in the light of the major approximation involving only 
nearest neighbour interactions. 


§ 4. StrREsS TO MOVE THE SUPER-DISLOCATION 
The stress to move the super-dislocation is calculated from eqn. 
(4), where (Brown and Herman 1956) 
Epp = l2NkTS,? et City vas (eis) 
and 
Ea V2ZNET IS a eee ee) 
where S, is the degree of order at X =a/2 (the slip plane). From eqn. 
(11) 
S,=Sosin [(P,-T)/T.P", | 
S2~ST,-T)L, for T>0sT,. | 
Substituting eqns. (13) and (14) in (4) 
T= NDS. AOs (6a oe ee LD) 
where 7 is the stress to move a super-dislocation. 7 is plotted against 7’ 
in fig. 4, for the case of CuZn where Sy was obtained from x-ray data by 
Chipman and Warren (1950), V =8-2 x 10! atoms/em?k= 1-4 x 10-8 ergs/°K 
and b=2-57x10-8cm. A very sharp maximum appears at 515°xK, 
where the yield point becomes 0-54 x 10®dynes/em?. In comparing the 
above theory with the experimental data it must be remembered that it 


(14) 
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was assumed that slip occurred in the absence of diffusion and also the 
effect of thermal fluctuations in reducing the yield point was neglected. 
These effects are expected to reduce the yield point but are not expected 


to cause much of a change in the temperature at which the maximum 
yield point is observed. 


Fig. 4 
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Theoretical curve of yield point versus temperature in region of long range 
order. 


§ 5. EXPERIMENTAL Data 


Two compositions of B-brass were used. of 48-6 and 49-9at.% Zn. The 
tensile specimens were 6-4mm in diameter and were both polycrystalline 
and single crystals. The yield points were determined on the Instron 
machine. For temperatures above room temperature the specimens were 
in air surrounded by a furnace, and below room temperature the specimens 
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were immersed in a coolant. No essential difference was obtained for 
single-crystal and polycrystalline specimens with respect to their 
temperature dependence of the yield point. Care was taken that only 
single phase material was used. The single phase 49-9at.% Zn was 
obtained below 300°C by heating to 750°c, water quenching, and then 
heating to 300°c for five hours. 
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Experimental curves of yield point and shear modulus of B-brass (48-6% Zn) 
single crystals versus temperature. 
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The yield point as a function of temperature is shown in fig. 5 for the 
48-6at.% Zn. Data by Ardley and Cottrell (1953) for B-brass of about 
the same composition are included. The variation of yield point with 
temperature may be divided into four regions. Region I is above T', and 
represents the behaviour of disordered B-brass : Region II (fig. 6) in the 


Fig. 6 
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Shows short range order effect above 7’,. 


immediate vicinity of 7', is probably related to the Sumino theory because 
Sumino predicted exactly the same kind of jump in yield point at 7',: 
Region III between 7’, and the minimum is associated with the present 
theory of the super-dislocation :.Region IV below the minimum is 
connected with the phase transformation found by Massalski and 
Barrett (1957). ; hoe 

Regions I and II are shown more clearly in fig. 6. It is important to 
note that 7 rises rapidly just before the temperature reaches 2,5 this rise 
is attributed to short-range order effect as presented by Sumino. It is 
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this short range order that also produces the tail in the specific heat curve 
(Dehlinger 1939). 

Ardley and Cottrell suggested that Region III was caused by the 
maximum in shear modulus (Good 1941) which is also plotted in fig. 5. 
If the shear modulus were responsible, then it means that a 6% change 
in @ can produce a 40% change in 7. This would not be possible if G 
acted by means of the internal stress where 7/G should be invariant as is 
clearly not the case. Another possibility is to assume that G effects the 
Peierls—Force in which case t might vary exponentially with G@ as indicated 
by Cottrell (1953). However, the slow change in order with temperature 
in the vicinity of 200°c does not indicate that G would behave anomalously 
in this region. It is suggested that the observed change in G is an effect 
of the yield point and not its cause, i.e. the change in G is a relaxation 
effect caused by dislocation movement. In fact, Clarebrough (1957) 
showed that a well-defined minimum in the internal friction existed at 
about 200°c. It is this relaxation in G@ [III] that would also account for 
the anomalous anisotropy observed in the elastic behaviour of £-brass 
(Good 1941). 

If one examines the phase diagram (Beck and Smith 1952) the 
composition of 48-6°% Zn appears to cross the two-phase boundary at 
about 200°C where the two-phase boundary at 200°C is obtained by 
extrapolation. In order to test the vital point of whether the maximum 
at 200° is associated with the Ba transformation, the 49-9 at. % Zn was 
tested. The composition of 49-9 at. °% Zn appears to cross the two-phase 
boundary at 75°c; if the maximum in 7 is associated with the B>a 
transformation, then Tmax should be shifted to a lower temperature for 
the higher Zn concentration. Figure 7 shows 7 versus 7’ for 49-9 at. % Zn 
where the maximum 7 is shifted to a higher temperature, 220°c, which is 
very close to the theoretical value of 242°c (fig. 4). In addition, the 
temperature of the minimum is shifted to a lower temperature, 50°C, as 
predicted by Massalski and Barrett, i.e. an increase in Zn content decreases 
the transformation temperature. 

For 48-6 and 49-9at.% Zn the maximum yield point is 0-33 and 
0-367 x 10° dynes/em? respectively as compared to the theoretical value of 
0-54 x 10° dynes/cm?. Since B-brass orders very rapidly, it is thought 
that diffusion effects and thermal activation are jointly responsible for 
the yield point being less than theoretical. 

Ardley (1955) has measured the yield point in AuCu, and found a 
maximum just below 7’,. This maximum seems to fit the Sumino theory 
better than the present theory. However, the maximum predicted from 
the present theory should occur ata higher temperature for AuCu, than 
for CuZn because AuCu, is more ordered close to 7',. Additional data 
showing a maximum in yield point versus temperature was obtained by 


eee 


} The experimental value was reduced by one-half because the specimen was 
polycrystalline. 
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Marcinkowski (1958) for Ni,;Mn. It is concluded that if yield point is 
measured as a function of the equilibrium degree of long-range order, a 
maximum will be observed. No data are known to exist that contradict 
the above conclusion for the cases where no change in crystal structure 
occurs during the long-range order transformation. 


Fig. 7 
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Yield point versus temperature of polycrystalline B-brass (49-99% Zn). 


§ 6. CONCLUSIONS 


The width of an anti-phase boundary was calculated. The stress to 
move a super-dislocation in the AB b.c.c. super-lattice was found to be 
NkT'S8,7/46 
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where S, is the degree of long-range order and 7’ is the temperature. For 
B-brass a maximum at 7'=515°K is predicted and the maximum is 
observed at 493°xK. The theoretical and experimental values of tmax are 
0:54 and 0:36 x 10° dynes/cm? respectively. 

A very sharp increase in yield point was observed at 7’, which agrees 
with the theory by Sumino. 

There is a minimum in yield point which shifted from 50 to 75°c when 
the composition was charged from 49-9 to 48-6 at. °/ Zn in accordance with 
the b.c.c. to close-packed phase transformation observed by Massalski 
and Barrett. 
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ABSTRACT 


Measurements are made in mixtures of aniline and cyclohexane at various 
concentrations and temperatures in the vicinity of the solution point of the 
velocity and attenuation of ultrasonic pulses. It appears that the large 
values of the attenuation observed in this region are due to the scattering 
of the radiation, which has been measured as to its angular distribution at 
various frequencies. 

Evidence from the extinction and directionality of sound signals in the 
mixture shows that the scattering is due to biphasal clusters which exist in 
the neighbourhood of the solution point and which are not set up by the 
passage of the sound itself. There remains some absorption which is of a 
molecular (viscous) nature. 


§ 1. INTRODUCTION 


MEASUREMENTS of the velocity and attenuation of sound in mixtures of 
liquids have been made fairly frequently since sources of ultrasonics became 
available and have shown in many cases a Maximum or minimum in both 
factors at a certain concentration. ‘This is so in liquid mixtures which 
form a complete solution. Our interest has been in those mixtures which 
in a certain range of concentration and temperature persist in two phases, 
whereas in another range they form a homogeneous solution. At certain 
concentrations the interface vanishes in a way similar to the disappearance 
of the liquid phase at the liquid—vapour critical point. By analogy, this 
phenomenon is known as ‘critical solution’ and the temperature as the 
‘critical solution temperature’. 

Interesting sound transmission phenomena have been noted in the 
neighbourhood of the critical state of vapours. It has been the purpose 
of this work to extend the study of these effects. Thermodynamically, 
the critical solution point of a pair of liquids is similar to the liquid-vapour 
critical point and provides an admirable opportunity of studying the 
one-to-two phase transition. In many cases the critical conditions occur 
at atmospheric pressure and a temperature near that of the laboratory. 

The present investigation has as its principal object a study (not hereto- 
fore attempted) of the mechanism of the ultrasonic attenuation in the 
Oe a ee 2 ee 
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critical region. One important question is, are the factors which produce 
the attenuation inherent in the mixture or are they set up by the passage 
of the radiation? Most of the previous investigators have made the latter 
assumption. 

There has been a tendency in this field to make a few measurements on 
a number of different mixtures rather than to study one in detail, as we 
decided to do. 

The absorption of sound in binary mixtures, is of course, closely related 
to the viscosity of the mixture. If there is no change of volume on mixing, 
the viscosity away from the critical state is linear with concentration, but 
if, as is often the case, there is a contraction, the viscosity shows a maximum 
at an intermediate concentration. This is so in alcohol—water solutions 
studied by Sette (1955). If, however, the concentration is kept constant 
while the temperature is varied and the liquids can pass from the double— 
phase to the single-phase type at a certain temperature, then the attenua- 
tion of sound rises to a value many times that which would be anticipated 
from the viscosity, as Chynoweth and Schneider (1951), also Alfrey and 
Schneider (1953), have shown, in aniline—hexane and triethylamine—water 
mixtures. Cevolani and Petralia (1952) found similar behaviour in aniline— 
hexane and phenol—water mixtures. When solution takes place, the 
velocity of sound has a value intermediate between that proper to the 
individual liquids. The experiments are usually carried out at ultrasonic 
frequencies from 1 to 30 Mc/s. 

In all these papers the data on attenuation and velocity were recorded 
but no experiments were carried out on the mechanism of the transition in 
relation to the ultrasonic propagation. 


§ 2. THe PRESENT EXPERIMENTS 


In the choice of a system for the study of critical solution, the following 
considerations were borne in mind: 


(a) The solution temperature must be neither too near nor too remote 
from room temperature. 


(6) The liquids must be obtainable in a good state of purity. 


(c) To minimize evaporation, the critical solution temperature must be 
well away from the boiling point of either liquid. 


The system aniline-cyclohexane seemed to fulfil these conditions best. 
The critical solution temperature is near 30°c, minimizing problems of 
thermal insulation and stability. The substances could be obtained pure 
in reasonable quantity. Cyclohexane boils at 83°c and aniline at 184°c, 
well above the critical solution temperature. 

Measurements of velocity and attenuation were made, using a pulse 
technique, at 1-5, 2-5 and 5 Mc/s. The composition of the mixture was 
varied and observations made for several degrees on either side of the 
critical solution point. 
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Following techniques which have been used for a number of years in this 
laboratory, a double transducer system of quartz—one to emit and one to 
receive—was employed. 

Pulses of 40 microseconds duration were derived from a gated oscillator. 
These were applied to the transmitting crystal which was carried on a 
movable piston. The electrical signal from the receiver, having an 
amplitude proportional to that of the incident sound, was, after amplifica- 
tion and detection, displayed ona cathode-ray oscillograph. As the crystal 
separation was reduced, the amplitude of the display was kept constant 
by the operation of a calibrated attenuator in the receiving circuit, its 
reading corresponding to the attenuation of the decreased path. 

A graph of the movement against attenuator reading gives the attenua- 
tion coefficient «’ in decibels per centimetre directly from the slope; a’ is 
related to the amplitude attenuation coefficient « by the relation 

a’ =20alogye or «’ =8-686a. 

Measurements of the velocity of propagation of the pulse, required in the 
calculations of scattering, vide infra, were made by mixing the signal at 
the receiving transducer with a comparable fraction of the (electrical) 
pulse applied to the sender. The magnitude of the resultant measures the 
relative phase of the components, which varies rhythmically with the 
separation of the transducers and so can be used to calculate the velocity. 
These experiments were made entirely in the Fresnel region of the radiation, 
i.e. for a separation of transducers less than R?/A where B is the radius of 
the quartz face (4 in.) and 2 the wavelength in the liquid. 

The vessel containing the liquids is a vertical glass tube 2 in. diameter 
and 6in. long. The lower quartz is pressed from outside against a thin 
aluminium disc, forming the base of the vessel; the upper is mobile and 
has the usual adjustments for setting it parallel to the other. The whole 
is then placed in a container 4} in. diameter and 8 in. long with heater 
coils and thermostat control. The latter is designed to hold the tempera- 
ture—read by thermocouple in the liquid—constant to 0-01°c over a 
period of 5 min. 


§ 3. PREPARATION OF LIQUIDS AND TESTS 


Benzene was used as test liquid. Adjustments of the mechanical part 
of the apparatus were made until a value for the attenuation in this liquid 
was obtained in agreement with the most precise of published results. 

Liquids of ‘Analar’ grade were used for the mixtures. As samples 
obtained on different dates from the suppliers gave consistent results, it 
was concluded that the purity of the liquids was sufficient to give reliable 
results. 

Each mixture was made up by volume, shaken in a bottle, then intro- 
duced into the apparatus, stirred to remove the bubbles, and the thermostat 
set in operation. A graph of micrometer screw reading against attenuator 
setting was plotted and the attenuation in decibels per centimetre calculated. 
Any graphs showing appreciable departure from linearity were discarded. 


222 
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The thermostat was set to allow a slight fall of temperature during the 
day. As the critical solution temperature was approached, the fall was 
slowed down by the emission of the latent heat of solution—a convenient 
effect since more measurements could be made in this important region, 
Measurements were made only at falling temperature. If the stratified 
liquids were left undisturbed, mixing did not take place until the temper- 
ature had risen a degree or more above the temperature at which separation 
took place, doubtless an influence of the size of the experimental vessel. 
Vigorous stirring was found to promote mixing, but often carried air 
bubbles under the edge of the upper crystal holder. It was therefore 
assumed that the temperature of separation obtained in a run at falling 
temperature was the true critical solution temperature and effects of 
hysteresis were not investigated. 

The experimental points obtained by a run of the above type usually 
fell on a smooth curve, indicating some degree of reliability. Intermediate 
points could be filled in on subsequent runs with quite good accuracy. It 
was estimated that the overall accuracy at higher attenuations ( ~ 2 dB/em) 
was 2%, at lower attenuations (>0-3 dB/em) 5%, the main uncertainty 
being the setting-up of the signal on the graticule of the oscilloscope. 

The concentrations studied ranged from 20 to 60% mols of aniline and 
the temperatures from 29° to 34°c by half degrees. 


Fig. 1 
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Velocity of sound in mixtures of various molar percentages of aniline in 
cyclohexane at 5 Mc/s. 


§ 4. VELOcITY RESULTS 


Measurements of velocity as a function of temperature for a range of 
concentration were made at 1-5, 2-5 and 5 Me/s. Figure 1 shows the 
results at the highest frequency, with the critical temperature indicated 
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as 30-7°c fora concentration of 33-39%. Below this temperature the mixture 
breaks into two layers, so that different results are obtained for the ultra- 
sonic wavelength according as the interferometer path lies entirely above 
or below the interface. 

Although the curves for the three frequencies did not exactly fall together, 
especially near the solution point, a fact which we ascribe to different 
mixing conditions for successive solutions, we have no definite indication 
of dispersion of velocity in this range of frequency. 


Total attenuation in mixtures at various frequencies. 


§ 5, ATTENUATION RESULTS 
While the ‘ triple point’ in the velocity always occurred at the solution 
temperature, the maximum of attenuation often occurred a little higher 
(up to 14°c higher). After analysis of many observations at four 
temperatures in this region, they were re -plotted as shown in fig. 2. They 
exhibit the rise of attenuation with frequency, when the mixture has a 
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concentration near 40%. (The pure liquids have values for the coefficient 
a ranging from 0-001 to 0-04 as the frequency rises in this range.) 

The previous workers in this field all ascribe this large absorption to 
relaxation. processes, mostly structural in nature, as though much of the 
energy in the sound beam were dissipated in forming and breaking up 
‘clusters’; indeed, some speak of ‘cluster relaxation’. Be this as it may, 
we have adopted the viewpoint that the clusters may be already present 
in the critical state—as they certainly are in carbon dioxide near its 
critical point, cf. Andrews’ well-known experiments—and may scatter the 
radiation out of the forthright beam, just as they scatter light in this region 
in the phenomenon of ‘ critical opalescence’. 

Scattering will occur wherever fluctuations in acoustical refractive index 
occur, having their origin in density variations or circulatory motions 
in the fluid. The former type is to be found in suspensions of solids or 
liquid globules in a matrix liquid (Allinson and Richardson (1958)); the 
latter in an atmosphere filled with eddies (Richardson (1950)). The 
scattering of sound from a rigid sphere, small compared to the wavelength 
A, was first calculated by Rayleigh and led him to the well-known law, 
which applies equally in optics. When the obstacle is large, the scattered 
energy is concentrated in a direction down the beam (behind the scatterer). 
At intermediate sizes (27a~ZA) the diffraction phenomena familiar in 
optics under the name of Mie make their appearance and withal the total 
energy scattered goes through sinuous variations (as A decreases) until it 
establishes itself at a value twice that intercepted by the scattering obstacle. 
Although the form of this curve is the same for all acoustic properties of 
the scattering material, the amplitude of these variations in the total 
scattered energy depends on the refractive index, to sound, of the material 
and on its density relative to the ambient medium. 


§ 6. ScaTTERING MEASUREMENTS 


In order to examine these aspects of the problem it seemed to us desirable 
to modify the apparatus so that the proportion of sound energy diffused 
out of the main beam could be measured. 

An ultrasonic beam is transmitted across a vessel containing the mixture. 
A receiving transducer or probe is mounted so as to allow it to be moved 
in a circle through whose plane the axis of the beam passes. Pulsing of the 
sound transmitter allows spurious echoes from the walls and the free surface 
to be distinguished by their associated time delay. The amplitude of the 
signal at any angle is read directly from its height on the oscilloscope 
display. 

A cross-section of the mechanical apparatus is given in fig. 3. The 
experimental vessel is of sheet steel, 6 in. deep and about 3 in. square in 
section. It is arranged to fit neatly into the thermostat cylinder whilst 
being suspended from a wooden base-board carrying the usual thermal 
insulation. A hole in the side of the vessel is covered with an aluminium 
diaphragm. Against the back of this the transmitting crystal is pressed 
by means of a leaf spring acting as the r.f. connection. 
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The receiving probe is also a quartz crystal carried in a brass holder 
having an aperture of lem. This size is chosen as a matter of compromise. 
In this work it is necessary that the angle at which the energy is scattered 
be known with precision. For a given wavelength, the receiver becomes 
more directional the larger its diameter, making a large-aperture receiver 
desirable. On the other hand, if the beam of scattered sound is of less 
width than the receiver, the latter will give an erroneous value of the 
scattered amplitude. The chosen diameter (1 cm) gives an angular half- 
width of 1-4° of the central maximum of its diffraction pattern at 5 Me/s. 
This diameter is less than half that of the transmitted beam and—it is 
assumed—of the scattered beams. The size of the chamber in these 
experiments is limited by the necessity to keep a close control on the 
temperature and by the high attenuation in the vicinity of the critical 
solution. 


Fig. 3 
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Consider a single scattering obstacle in a plane sound beam. Let the 
intensity of scattered sound at angle 6 to the beam (as measured by a probe) 
be I». Assuming the intensity of scattered sound is symmetrical about 
the beam axis, the scattered energy dH at angle @ is given by d# =I[(gdQ 
where dQ is the element of solid angle between 6 and 6+d0. The total 


energy 
0 


I, being the central maximum and f(@) the distribution with angle. 
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Now let the experiment be repeated in a pure liquid having the same 
velocity of sound so that the pattern due to diffraction alone will be 
unchanged. The total energy 


E,=2nly | " #7(6) sin 6 d8, 
0 


I, being the intensity in the central maximum. 

Both the central and azimuthal beams are subject to molecular dissipa- 
tion which may not be the same for the mixture and the comparison liquid 
but does not affect the polar curve in the first order. (The effect of attenua- 
tion on the theory of diffraction has been studied by Alblas (1957).) If 
therefore the central maximum is adjusted to the same value in each case, 
the energy scattered between the crystal and a plane of scansion at a 
distance x can be expressed by a coefficient 


1 is 
ot od 9g 8G. 


where P,, and Q, are the values of the integrals 
[ro nodosa [ro sin 6 d6 


at the distance 2 There is of course an implicit assumption that the 
pure liquid does not scatter. This might be tested by plotting the sound 
field round the transducer in a gas like argon for which the absorption is 
given by the classical formula of Kirchhoff and adjusting the pattern for 
the difference in sound velocity. There is also the ever-present difficulty 
that the scattering is multiple whereas the theory assumes single scattering. 
(Multiple scattering has been studied by Fejer (1953) but his theoretical 
results for a <A are not applicable to this case.) 

If the transmitter behaved as a piston source it might be thought that 
interference between the contributions at the receiver from all the scatterers 
within the volume which the probe ‘views’ would occur, but in practice 
the width of the probe seems to obliterate such interference, except at the 
highest frequency. The scattered distribution could therefore be mapped 
out directly without undue labour. In the 5 Mc/s case the pattern con- 
sisted of many peaks whose precise delineation would have been laborious 
and unnecessary. The following procedure was therefore adopted. The 
amplitude at the top of each peak was noted. A plot of these readings 
against angle was similar in form to one in which the peaks and valleys 
were smoothed out—i.e. in which interference was allowed for—except 
that compared to other curves the vertical scale was doubled. The 
value of # evaluated from a graph of J(, sin 6 by mechanical integration 
was therefore divided by 4. This procedure however obliterated the 
central peak but a value obtained from the experiments on non-scattering 
liquids could be inserted to compensate for this, if required. 
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§ 7. ScarrEeRING REsutts 


The angular distributions of sound were obtained and the scattering 
attenuation estimated at frequencies of 5, 2-5 and 1:5 Mc/s in mixtures 
containing 22-5, 33-3, 43-8 and 53-9 molar per cent of aniline. The measure- 
ments at 5 Mc/s occasioned little difficulty apart from the complicated 
nature of the distribution pattern, presumably due to interference. At 
2-5 Mc/s the distribution took on more the aspect of a broadened central 
beam, the angular spread not being more than 4° or 5°. At 1-5 Me/s this 
spread was greatest, although it was still possible to distinguish between 
scattering and non-scattering liquids with certainty. 


Fig. 4 


Scattered attenuation at various frequencies. 


The overall scattering measurements are presented in a three-dimen- 
sional plot (fig. 4) of the same type as that in fig. 2. 

At 22-5% aniline concentration the scattering attenuation almost 
keeps pace with the total sound loss throughout the critical region. At 
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Temper- 
ature °C 


22-56% 

aniline a Ne 
0-15 | 0-11 
0-15 | 0-11 
0-155 | 0-11 
OFT WW OsZ 
0-185 | 0-14 
0-26 | 0-145 


1281 
1284 
1288 
1294 
1296 


aniline 
34 
33 
32 
31 
30:5 
30 


33°3%, however there is a falling-off in scattering at a slightly higher 
temperature than the attenuation peak, the scattering attenuation not 
recovering its former value until around the solution point. With the 
43-8% mixture, the behaviour is most striking. The scattering attenua- 
tion, having kept pace with the total loss, falls sharply some 2°c above the 
solution temperature reaching a very low minimum at a temperature 
almost identical with that of maximum total attenuation. There are 
signs of an increase in scattering at temperatures below the solution 
temperature, but it is difficult to decide in which of the two liquid layers 
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the measurements have been taken and such results may be misleading. 
The 53-9% mixture displays little scattering, but a peak is noted near the 
temperature of maximum attenuation together with an increase through 


Fig. 5 
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Angular distribution of scattered sound at 5 Me/s. 


the solution temperature. One of the matters of note about the curves 
is the way in which the scattering and total attenuation keep pace at the 
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high temperature end of the range. When the scattering coefficient is 
subtracted from the total coefficient in this region the remainder is of the 
order of that found in the separate liquids, viz. aniline o/f?=50, cyclo- 
hexane a/f%=200, x 10-1’. 

The individual results are shown, as total attenuation and total scattering 
intablel. Although both reach a maximum near the solution temperature, 
the two peaks do not always fall together. 

There are several methods of analysing the scattering curves of a mono- 
disperse suspension of spheres to derive the particle size in terms of the 
wavelength, the relevant parameter being y= 27a/A: 


(1) From the variation of scattered intensity with angle using tables 
such as those compiled by Lowan (1946). 


(2) From the variation of extinction in the direct beam, supposed to be 
deprived of energy by scattering, with wavelength or frequency. 


Curves of the relative intensity of sound at various angles from the 
beam axis for certain temperatures and compositions are given in fig. 5. 
The relative intensities plotted are the average values for the required 
angle on both sides of the central maximum. In a medium without 
scattering the diffraction pattern of the probe falls to zero (first minimum) 
at an angle of 1-5° at 5 Mc/s. Viewing our curves in the light of this fact 
demonstrates beyond any reasonable doubt the existence of the scattering 
phenomenon. 

Maxima and minima in the polar curves are noted, particularly 
when the mixture contains one third of aniline (molar), indicating 
some degree of uniformity of size of the clusters. From the position 
of these maxima and minima and from the polar curves in general, it 
should be possible to deduce the mean cluster size, if the refractive index 
m is known. The velocities of sound in the pure liquids are 1250 and 
1660 m/s respectively giving m=1-33—as for rain drops in air to light 
waves. The calculation of the angular distribution can be made from the 
Rayleigh—Gans theory if the refractive index is sufficiently near to unity 
for the scattered field to be solely due to the phase difference between the 
wave which has passed through the ‘sound-soft’ sphere and that through 
the external medium when they reach the rear of the obstacle. Otherwise 
the calculation of the field is difficult but Hart and Montroll (1951) have 
been able to extend the Rayleigh—Gans calculations to refractive indices 
in the range m?—1<y—1 which would cover all the cases encountered in 
these mixtures. In what follows it appears that this cluster size is com- 
parable with or, in some cases, greater than the wavelength. When y is 
considerably greater than unity the distribution with angle takes the form 


I,0cl+ cot? 0/2J,?(y sin @) 


but in the awkward region where y=1, the tables constructed by Lowan 
for the corresponding light case must be used. 
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Pursuing this course for some trial values of y we find, for instance, that 
for y=3, the first maximum should occur at about 15°, which corresponds 
to the distribution at 30-2°c in the mixture of one-third aniline, while 
the steeper fall in 53-99% aniline at a nearby temperature would indicate 
y=5, with a much smaller value at 36°c, 


Fig. 6 


22-5 My ANILINE 33-3 %o ANILINE 


53.9 fp ANILINE 


6 6.2 6.4 6:6 6 6.2 6.4 6.6 
log f 


Frequency dependence of scattered sound. 


A log-log plot of scattering attenuation against frequency for various 
temperatures and concentrations is given in fig. 6. For the curves away 
from the absorption region the slope is round about 2, indicating a scattering 
proportional inversely to the square of the wavelength. 
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In the case of the 53-9°, concentration, the power changes in sign at a 
higher frequency, while in the 33-3% case it increases with decreasing 
frequency, approaching the value 4 appropriate to scattering by obstruc- 
tions much smaller in size than the wavelength. Such behaviour of the 
extinction coefficient—fraction cut-off the direct beam—recalls the 
maxima and minima of transmission as A is varied and whose position is 
again a function of y. It is usual to express the transmission in terms of a 
scattering cross-section (per unit area) A as it varies with y(=27af /c, with ec 
the velocity of sound, in this case). A graph of this parameter for a re- 
fractive index m= 1-33 has been calculated by Hart and Montroll (1951) 
and is shown in fig. 7 for comparison with fig. 6. 


Pigs 


0.5 1-O 1.5 
log 
Variation of extinction due to scattering with parameter 27a/A for refractive 
index, 1-33. 


We have made the comparison by choosing for each temperature a 
value of y to which each of the four sets of curves in fig. 6 of «, against f can 
be best matched in slope and turnover. Knowing the frequency and 
velocity of sound we then calculate from this value of y the ‘ cluster radius’ a. 

Here are some typical results at 2-5 Mc/s and 31°. 


Table 2 


Mol % aniline] + estimated A (mm) a (microns) 


9) 0-50 120 


0-51 150 
0-53 250 
0-54 400 
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It may well be that the value we have assumed for m is too high. There 
is no reason to suppose that the globules should be formed of pure aniline 
in a matrix of pure cyclohexane. A smaller refractive index would shift 
the first peak of the curve in fig. 7 in the direction of higher values of y 
in the ratio A/(m—1), but in the absence of knowledge of the coefficient 
of mutual diffusion of the two liquids, we cannot quote a precise value for m. 


§ 8. SCATTERING PROCESSES 


We should hesitate to aver that, within the limits of accuracy of the 
scattering measurements, the whole of the reduction in intensity in the 
direct beam is to be ascribed to scattering. After subtracting that which 
can be ascribed to shear viscosity, there remains an amount, of the same 
magnitude as that in the pure liquids, which must be put down to relaxation 
or bulk viscosity. 

We now revert to the question we raised earlier; do these clusters exist 
already in the liquids when they approach the critical solution temperature 
or are they caused by the passage of the ultrasonics? Ernst (1947) in his 
investigations on the velocity of sound—measured by optical diffraction— 
in another binary mixture mentions the appearance of a ‘brick wall 
structure’. He shows a photograph of this appearance but he does not 
give a scale, though it appears that the cells have a diameter of the order 
of the ultrasonic wavelength and they seem to us to correspond to the cells 
of Bénard which can be set up in an unstable liquid column as small scale 
copies of the ‘mackerel sky’, a cloud formation visible under similar 
conditions in the atmosphere. He remarks that the formation once set up 
could persist for hours. We have not observed this honeycomb structure 
in our mixtures but then we are putting in the energy only weakly 
and occasionally, Ernst used continuous waves of sufficient intensity to 
diffract light at the wave front. 

Though it is not possible to see clusters at the solution point with.the 
unaided eye because the relative optical refractive index of the liquids 
is only 1-11, we have successfully applied the schlieren technique to this 
problem. The schlieren apparatus was set up for the observation of 
phenomena due to convection in a large tank of water. For the present 
purpose a small glass-walled cell about 3-4 cm square and 7 cm deep con- 
taining the mixture was set up vertically in a large tank of water (acting as 
a thermostat) and in the way of the light beam of a schlieren apparatus. 
Under the schlieren view and in the absence of any sound propagation 
the clusters were clearly visible moving about, though they had rather 
the appearance of striations than globules, at all temperatures from 30° to 
33°o, and recalled the photographs made by Graham (1934) of a layer of 
liquid heated from below and sheared. When a photograph was taken 
of our striations it was possible to estimate their size e.g. at 32-5° there 
were 25 per cm length, at 33°0, 30 per cm (cf. fig. 8), Pl. 79. 

These inhomogeneities are the evident cause of the scattering which we 
have encountered in the ultrasonic experiments. Itis intended to continue 
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this work on sound scattering by density inhomogeneities to the class of 
fluids investigated by Graham, i.e. to set up temperature gradients in a 
single fluid and observe the effect on an ultrasonic beam transmitted through 


it. 
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ABSTRACT 


A technique is described whereby the surface of metal specimens can be 
examined by transmission electron microscopy. The original specimens are 
thinned by electropolishing from one side ; thus the surface contours, e.g. 
slip steps, are preserved on one of the surfaces. It is shown that dislocations 
are revealed associated with the surface slip lines; the step heights and 
surface contours such as those of fatigue striations are observed at a high 
magnification due to a taper section effect. 


§ 1. INTRODUCTION 


In recent studies of the dislocation arrangements in deformed metals by 
the transmission electron microscopy technique (Whelan ef al. 1957, 
Hirsch et al. 1958) bulk specimens were deformed and subsequently 
thinned by electrolytic polishing from both sides of the sheet. In this way 
the dislocation arrangement in the interior of the metal was revealed. Itis 
of particular interest however to study the dislocation arrangement and 
other defects near the surface of the specimen. Accordingly a technique 
has been developed whereby the specimens are prepared by thinning from 
one side only, so that one original surface is retained. In this way surface 
markings, such as slip lines, can be seen at the same time as the dislocations 
and other defects just below the surface of the metal. The experiments 
have been performed on stainless steel with the following objects in view: 


(1) To relate the slip lines to the dislocations just below the surface of the 


metal. 

(2) To examine whether the dislocation arrangement near the surface is 
similar to that in the interior. 

(3) To study the dislocations and other defects and to search for possible 
cracks at the surface of specimens deformed in fatigue, where the surface 


is known to be particularly important. 


§ 2. EXPERIMENTAL TECHNIQUE 


Specimens of annealed stainless steel with the same composition as that 
used in previous investigations (Whelan et al. 1957, Whelan 1959) and 
some specimens of a similar austenitic stainless steel (18/8 type) were 
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electropolished in a bath of 40° sulphuric acid and 60% phosphoric acid 
and then either deformed 10% in tension or fatigued for 107 cycles in reverse 
plane bending with zero mean stress. The thicknesses of the tensile and 
fatigue specimens were 0-:005in. and 0-012in. respectively. The fatigue 
machine was a commercial Avery plane bending machine, this mode of 
deformation being chosen because it is comparatively easy to maintain zero 
mean stress in a thin specimen. The fatigue stress was chosen so that 
10’ cycles corresponded approximately to 10% of the life of the specimens. 

After deformation the specimens were coated on one side with ‘ Lacomit’, 
a commercial acid resistant varnish, and thin films were then prepared by 
electropolishing from the unprotected side towards the coated surface. 
When the test pieces were thought to be sufficiently thin, specimens were 
cut off in the usual way (Tomlinson 1958) and the varnish dissolved with 
‘Lacomit Remover’. The thin films were then examined in a Siemens 
and Halske Elmiskop I operating at 80 or 100 kv. 


§ 3. EXPERIMENTAL RESULTS 
3.1. Observations of Slip Lines 


The surface slip lines can be observed on the electron micrographs in two 
distinct ways: 


(1) They appear as lines along the directions of the intersections of the 
slip planes with the previously protected surface, in much the same way as 
slip lines on surface replicas (e.g. Wilsdorf and Fourie 1956, Mader 1957). 


(2) Theslip lines are associated at the edge of the thin foil with steps which 
correspond to the slip line steps magnified by a taper section effect. 


3.1.1. Contrast effects 


Figure 1 (Pl. 80) is a micrograph of a stainless steel specimen deformed 
in tension showing a set of prominent slip lines crossed by a set of fainter 
lines corresponding to slip on another system. The transmitted intensity 
is the same on either side of a slip line, but the line itself shows very pro- 
nounced contrast. Sometimes the contrast of the lines is darker than the 
background, and sometimes lighter as for example in fig. 2 (Pl. 81), which 
is a micrograph from a similar specimen. Some of the lines on fig. 1 show 
both dark and faint white contrast. All these contrast effects are very 
similar to those observed on slip traces produced by movement of disloca- 
tions in thin foils (Hirsch e¢ al. 1956, Whelan e¢ al. 1957). By analogy 
with the interpretation of the slip traces, the contrast of the slip lines in 
the present case is considered to be due to strains in the oxide film on the 
specimen surface. However, whereas slip traces can be formed as the result 
of the movement of single dislocations, many of the slip lines in figs. 1 and 
2 are due to the passage of a number of dislocations. Another important 
difference between slip lines and slip traces is the fact that slip lines only 
occur on one side of the specimen, while slip traces consist of pairs of lines 
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corresponding to the top and bottom surfaces of the foil. This can be 
seen very clearly in fig. 3 (Pl. 82), where dislocations on the slip planes are 
associated with one line corresponding to the protected surface. Other 
clear examples can be seen at A in figs. 1, 2. 

The question arises whether the slip lines are the result of dislocation 
movement during deformation, or whether they reflect a rearrangement of 
the dislocations in the thin foil (thinned specimen). Movement of dislo- 
cations after thinning would result in the formation of slip traces (i.e. 
pairs of lines). These are observed only very rarely. In some areas, 
for example at B in fig. 1, the dislocations (below B) appear to be associated 
with slip lines at both ends, as if the two slip lines constituted a slip trace. 
Close examination however reveals that the dislocations are associated 
exactly only with one line; thus the set of dislocations below B belong to 
the slip line on their left, while the set of dislocations to the right of B to 
the other slip line. Furthermore the asymmetry of the contrast of the 
slip lines in fig. 1 is always in the same sense (i.e. white line on the right of 
black line), whereas for the pair of lines constituting a slip trace the 
asymmetry would be in the opposite sense (for examples see Whelan et al. 
1957). We conclude therefore that almost all of the lines correspond to 
slip lines produced either during deformation or during polishing. (The 
same is also true for specimens thinned from both sides where slip traces 
are hardly ever observed immediately after insertion in the microscope.) 
The strain at the oxide metal interface must therefore persist over periods 
of hours or days. This is quite consistent with observations on slip 
traces in stainless steel; examination of slip traces in the same specimen 
after an interval of several days has shown that the traces in fact persist 
over this period of time. 

It is possible however that some rearrangement of dislocations takes 
place during polishing. It seems unlikely that much rearrangement 
could occur, partly because the observed networks are often very 
complex, and partly because experiments on the movement of dislocations 
have shown that the dislocations are strongly pinned at the surface; 
thus while during thinning some of the obstacles blocking the dislocations 
might be released, the surface pinning force becomes more important, and 
furthermore, the range of elastic interaction between dislocations decreases 
(because they become shorter); both these effects would tend to prevent 
rearrangement of the dislocations. Nevertheless, the possibility of some 
rearrangement cannot be ruled out; for example, in metals in. which 
cross-slip is possible, when the foil is thin those dislocations in pile-ups 
which cross the foil from top to bottom could change into the screw 
configuration and cross-slip out of the slip plane, thus dispersing the pile- 
ups. This is clearly not the case in stainless steel. It is nevertheless 
planned to test for possible rearrangments of dislocations during thinning 
in this and other metals by close correlation between surface slip lines and 
dislocations inside the specimen, and by comparing surface replicas before 
and after thinning of specimens. 


342 
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In some areas slip lines can also be recognised as a result of a different 
contrast mechanism. In fig. 4 (Pl. 83) the transmitted intensity changes 
abruptly on either side of a slip line, and this is due to a change either in 
thickness of the specimen produced by the slip step, or in orientation of 
the crystal due to a set of dislocations on the slip planes. It should be 
noted that the transmitted intensity is a rather complicated (oscillating) 
function of thickness and of orientation; nevertheless abrupt changes of 
contrast as in fig. 4 can be taken as an indication of changes of either of 
these parameters. 

The same effect may be seen in a different way in fig. 3, where a thickness 
contour runs across the slip lines. At each slip step the contour shifts 
(for example at A, B) as would be expected from a change of thickness, or 
of orientation. It should be noted that the former interpretation is valid 
in regions free of strain; on the other hand a set of screw dislocations on a 
slip plane would tend to twist the parts of the crystal on either side of the 
slip plane relative to each other, and this change of orientation would 
lead to a shift of the extinction contours. In fig. 3 however there appear 
to be rather few dislocations in the neighbourhood of the contour 
suggesting that the shifts are due to thickness changes. 


3.1.2. Steps at the edge of the thin foil 


In fig. 3 slip lines can be seen to be associated with steps at the edge of 
the electron microscope specimen. ‘This edge is the line of intersection 
of the polished surface and the protected original surface which make small 
angles with each other at the edge. The edge of the thin foil is clearly wedge 
shaped. Because of the nature of this geometry a change in thickness of the 
specimen at a slip step gives rise to a much larger step at the edge of the 
thin foil. This magnification is due to the well-known ‘taper section’ 
geometry applicable in this case. It should be noted therefore that the 
steps at the edge of the thin foil correspond to a (considerably magnified) 
increase in thickness of the specimen, i.e. they correspond to the component 
normal to the foil of the slip step. Figure 5 shows the geometry of this 
effect. The plane ABDF is considered to be normal to the direction of the 
electron beam. BODE is the slip plane, the projection of BE, i.e. BF is 
the slip line and BC is the step at the edge of the foil. It is often possible 
to determine directly from the micrographs the taper angle in the active 
slip planes; thus from fig. 5 it follows that 1=a¢é (for small ¢) where 
Lis the slip step height measured in the slip plane and normal to BE, a is 
the length BC of the step at the edge of the foil measured along the inter- 
section of the slip plane with the polished surface (ie. BD), and ¢ is the 
taper angle (assumed small) in the slip plane. Further, /=aysec 0, 
where 7 is the projection of ¢ on to the lower surface of the foil and @ 
is the angle between the slip plane and the plane of the foil. It should be 
noted that the foil could also be mounted in such a way that the plane 
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ABE is normal to the direction of the electron beam. Since ¢ and are 
small in practice the results of the above analysis apply. 

The taper angle ys can be measured in a number of ways. For example 
in fig. 3 the dislocations become progressively shorter towards the edge, 
giving the taper angle ys. Or ys may be determined from a suitable slip 
trace, for example one nucleated at the edge (see Whelan et al. 1957); 
again, the variation in taper can sometimes be determined by measuring 
the thickness at various distances from the edge, either from the width of 
slip traces or lengths of dislocations in known slip planes. 

In fig. 3 the shifts of the extinction contour at A and B correspond to 
the steps at the edge; this confirms the view expressed above that in the 
absence of strain the shift of the contours is due to a change in thickness. 
Since the shift of the contours is approximately equal to the length of the 
steps at the edge, the taper is probably reasonably uniform between the 
contour and the edge. 


Fig. 5 


A B 


Geometry of taper section showing slip step magnification produced 
at the edge of the thin foil. 


With an instrumental magnification of 100 000 times and a taper magni- 
fication of 10 times a millimetre edge step would correspond to 104. Thus 
the technique is obviously a powerful one for studying the distribution 
of slip. Both the slip step height and the sense of the slip on different 
slip planes can readily be determined. The extreme edge of the foil 
would be expected to be rounded off during polishing. The only effect of 
this will be for the edge steps also to be rounded off. In fig. 5, ABC is the 
idealized shape. Rounding off will occur at B since polishing takes place 
along AB and BC and this effect can be seen in fig. 3. Nevertheless the 
length a can always be measured by extending the straight parts of the 
edge to intersect DC produced. 
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3.2. Correlation of Number of Dislocations on Slip Plane and 
Slip Step Height 

In fig. 3 the number of dislocations on the slip planes can be counted and. 
knowing the taper angle the offset expected from these and that measured 
by the taper section technique can be compared. In the case of plane B, 
where the dislocations can most readily be counted we have 12 dislocations 
and hence 1=304 for the Burgers vector nearest the foil normal. For 
this step ax~2-5mm at a magnification of 70000 times, 0=72° and 
~2-0°; hence 1~40+204. These two estimates of / are consistent 
with each other, but considerably better accuracy is required to make the 
correlation precise. 

Figure 6 (Pl. 81) suggests that it may be possible to make a similar 
comparison in another way. At A the slip steps on one system are offset 
by slip on a second system. Now between the two offsets indicated there 
are a number of dislocations on the crossing system. The difference in 
the offsets though small is measurable and can be compared with the 
displacement expected from the ten dislocations between them. The 
difference between the offsets is about 40 A which is twice as great as the 
value corresponding to the observed dislocations. This difference is 
presumably due to slip occurring on one of the offset systems after the 
other and after the crossing system had begun to operate. 

It is possible to check the accuracy of the taper section measurements by 
comparing them with the offsets on slip lines where slip has occurred on 
two systems. In fig. 4 offsets and edge steps can be compared but, because 
there is a choice of Burgers vectors, no definite conclusions can be drawn, 
although the two types of measurement are found to be consistent with 
each other. 

The results described in this section are of a preliminary nature, and are 
intended to show only that it is possible with this technique to make 
quantitative estimates of the magnitude of features on the atomic scale. 


3.3. Dislocation Arrangement at the Surface after Deformation 
in Tension 

The arrangement of the dislocations at the surface is found to be very 
similar to that observed previously in the interior of the specimen (Whelan 
et al. 1957). Thus the dislocations are generally arranged in pile-ups or 
networks on the slip planes (figs. 1, 2). These networks are the result of 
interaction of the piled up dislocations with others on intersecting systems. 
Figure 2 shows many extended and contracted nodes of the type discussed 
in detail by Whelan (1959). It appears that some of the dislocations taking 
part in the network A are on a slip plane slightly displaced from that 
containing most of the dislocations. Although no systematic statistical 
average has been obtained from the micrographs, the networks or pile-ups 
often contain about 20-30 dislocations corresponding to slip displace- 
ments ~ 50-754 in agreement with the order of magnitude of the size 
estimated from the edge steps (see figs. 3, 4). 
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Examination of the sense of the edge steps shows at once that in general 
this is the same for most of the slip steps after unidirectional deformation 
(see figs. 3, 4). Taking the average displacement to be ~504, and the 
spacing of the slip lines ~ 10004, the average shear strain is found to be 
~5%, i.e. of the correct order of magnitude (the deformation of the 
specimens corresponding to figs. 3, 4 was 10%). 

While many of the dislocations are arranged in networks at the surface, 
as for example in fig. 2, it should be noted that if an etch pit technique had 
been used to reveal the dislocations, a sequence of pits resembling that of 
an irregularly spaced pile-up would have been observed. It is clear 
therefore that the interpretation of an array of etch pits in terms of a 
detailed arrangement of dislocations has to be carried out with caution. 


3.4, Some Observations on Fatigued Specimens 


The surface of a metal is particularly important in the study of fatigue. 
The first fatigue cracks are often if not always formed at the surface. 
In a number of face-centred cubic metals these cracks have been shown 
to occur during the first few percent of the life in the broad slip bands 
known as striations which develop during fatigue (for review see Thompson 
and Wadsworth 1958). Itis ofinterest to study with the present technique 
the distribution of slip and the nature of the defects in the striations, the 
surface contour across them and, if possible, to observe the cracks at an 
early stage of their formation. 

Figure 6 shows the appearance of the surface of a fatigued specimen. 
Area B shows slip lines and a dislocation arrangement not unlike that in 
a unidirectionally deformed specimen. At C however, there is a striation 
consisting of many closely spaced slip lines. There is also a very high 
density of dislocations in this region ; the resolution is however too poor to 
make any estimate of this. Another feature often observed is the presence 
of patches of stacking faults in the slip bands, as at D, or in the striations. 

Figure 7 (Pl. 84) shows the distribution of slip near the edge of the 
thin foil. In contrast to the micrographs obtained from unidirectionally 
deformed specimens the sense of the slip steps alternates frequently. 
Another distinguishing feature is the fact that some quite distinct slip lines 
have no steps associated with them, for example at A. While this could 
be due to slip with a Burgers vector in the plane of the foil, it is likely that 
the Burgers vector corresponding to the slip lines in the same area will be 
the same. We suggest therefore that these lines are the result of to and 
fro slip either on the same plane or on closely neighbouring planes. Some 
of the lines also show lighter contrast (for example at B) of a type not 
observed in unidirectionally deformed specimens. While the inter- 
pretation of these lines is not yet clear it is possible that they correspond 
to the early stages of cracks along the slip lines. 

Figure 8 (Pl. 85) shows a very pronounced striation. Here the surface 
contour of the striation is revealed by the steps at the edge of the thin 
foil. The taper magnification is approximately 10 times, the magnification 
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of the micrographs is 45000 times, and hence the magnification of the 
steps is approximately 450000 times. This is to be compared with a 
magnification of 20 x 1000 times achieved in recent taper section experi- 
ments of Wood and Segall (1957) and Wood (1958), who used an optical 
method. Itisclearfrom this micrograph that slip can take place in opposite 
directions at spacings as close as ~504. Both projections and fissures 
can be observed. Some of the former project above the mean surface 
level (i.e. in a direction normal to the plane of the foil) by about 3004. 
(It should be remembered when examining this micrograph that the 
magnification at the edge in the direction of the arrow is about ten times 
that at right angles to it. Hence the height of the projections is about 
equal to the width of their base.) The peaks undoubtedly represent genuine 
projections above the surface and may correspond to the early stages of the 
extrusions described by Forsyth and Stubbington (Forsyth and Stub- 
bington 1955, Stubbington and Forsyth 1957). While the fissures 
correspond probably to intrusions, it is possible that they might be due to 
preferential polishing at a slip plane. The same argument applies to the 
white lines which can be observed and which sometimes end in fissures. 
In principle it might be possible to determine whether the white lines 
correspond to depressions in the protected surface or polishing groves on 
the other surface; in the former case they should appear as single lines; 
in the latter they should be paired with a slip line. There is little doubt 
that some of the white lines in fig. 7 represent actual fissures since they 
occur singly. In fig. 8 the slip lines are so closely spaced that it is difficult 
to be sure of the interpretation. Whatever the origin of the white lines 
or fissures might be, there is no doubt that they correspond to structural 
features not present in unidirectionally deformed specimens. Close 
examination of the white lines in fig. 8 also reveals that they consist of a 
series of closely spaced short white lines running transversely across them ; 
the origin of this structure is not understood. 


§ 4. CONCLUSIONS 

The preliminary results obtained so far with the technique described 
in this paper for studying the surface of a metal, show that it combines 
the advantages of the transmission and replica methods of electron 
microscopy. It has been shown that in stainless steel it is possible to see 
the dislocations in the metal associated with the slip lines. The distribu- 
tion of dislocations at the surface is similar to that in the interior of the 
metal. The height of the slip steps and surface contours, for example of 
fatigue striations, are revealed at a high magnification by the steps at the 
edge of the foil due to a taper section effect. It appears that the technique 
is likely to prove particularly useful for correlation between dislocation 
arrays and slip lines, and for the study of the distribution of slip, fissures, 
projections, and dislocation arrangements in fatigued specimens, 
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ABSTRACT 


The structures of some fifty alloys of silver, manganese, and aluminium 
have been studied with a view to explaining the magnetic properties of these 
alloys. Approximate equilibrium diagrams were made at 400°c and 650°C, 
and the alloys were also examined after quenching from the latter temperature. 
Tt follows from this examination that ferromagnetism in alloys in the vicinity 
of the atomic composition Ag,;MnAl is not associated with any equilibrium 
product, but arises from a finely dispersed nucleation structure, which grows 
during ageing into a Widmanstitten pattern and eventually transforms above 
the Curie point into B-manganese, rendering the alloy non-ferromagnetic. 
The ferromagnetic precipitate was extracted and identified by electron 
diffraction as having an ordered hexagonal structure with a=—2-94 4, 
c/a=1-61. This result is used to discuss the variation of magnetization with 
temperature in these alloys, and also the origin of their very high coercivities. 
A brief discussion of the valency of manganese in these structures is also 
included. 


§ 1. INTRODUCTION 


Tue silver-manganese—aluminium alloys are ferromagnetic and belong 
to a series known as the Heusler alloys, after their discoverer (F. Heusler 
1903). The original compounds, discovered in 1898, were copper-based 
alloys, and it was later found (O. Heusler 1928) that the maximum intensity 
of magnetization arose in alloys, centred around the composition Cu,MnAl, 
which had been quenched from a temperature well above the Curie point, 
and then subjected to a low temperature ageing treatment. The structure 
associated with the ferromagnetic phase was not solved until Bradley and 
Rogers (1934) concluded that it was an ordered structure based on a body- 
centred cubic cell, the atoms in the body centre being alternately manganese 
and aluminium, and the corner atoms copper. 

The solution of this structure presented a fundamental difficulty, because 
the x-ray scattering factors of copper and manganese are virtually identical. 
Prior to this, Potter (1931) had endeavoured to overcome this by sub- 
stituting monovalent silver for monovalent copper, but in doing so, he 
found that the maximum intensity of magnetization arose not at Ag,MnAl, 
butat Ag;MnAl. Figure 1 shows Potter’s results, redrawn to give contours 
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of equal saturation intensities per cubic centimetre. The cross marks the 
composition Ag,MnAl. 

Potter discovered three further facts. The maximum intensity of mag - 
netization was achieved by chill-casting and ageing the alloy for several 
days at 250°c. Secondly, the magnetization was rapidly destroyed by 
heat treatment above the Curie point (~350°c), and finally the x-ray 
powder patterns taken from this alloy in its magnetic state showed only the 
pattern ofa silver solid solution. Thus, no evidence could be obtained for a 
phase isomorphous with the structure of Bradley and Rogers. 


Fig. 1 


40 


ite) 20 30 
Mn {at %o.) 


Saturation intensities /cc. of Ag-Mn-Al alloys (atter Potter ) 


Magnetization contours in e.m.u./cm® of the ternary alloys (after Potter (1931)). 


The alloy Ag;MnAl, which has been given the name ‘Silmanal’ is charac- 
terized by a low intensity of magnetization, but has a very large coercive 
force. Values of 50000e are common, and more than 80000e can be 
reached under certain conditions. Apart from advances in the technical 
applications of these alloys, however, no further fundamental investiga- 
tions were made on the structures of the alloys, with the exception of a brief 
report by Geisler (1951) which will be discussed in §6. A detailed survey 
was therefore undertaken in this Department. 
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§ 2. MarERIALS AND APPARATUS — 


The alloys used in this investigation were prepared from fine silver grain 
(99-99%, Ag), Johnson-Matthey H. S. manganese, and high-purity alumi- 
nium. The materials were melted together in an argon-are furnace, and 
homogenized in evacuated silica tubes at a suitably high temperature close 
to the solidus of the alloy concerned. For the purposes of identification in 
the course of this paper, the alloys will be enumerated in terms of their 
atomic compositions, although this does not imply the presence of any 
intermetallic compounds, nor that the actual composition as shown by 
subsequent analysis is that indicated by the formula. 

Microscopic examination was made on specimens polished either with 
fine alumina or liquid silver polish (‘Silvo’), and etched in most cases in a 
5° solution of potassium cyanide and ammonium persulphate, although, for 
alloys low in silver, an acid ferric chloride etch appeared more satisfactory. 
The preparation of these alloys for photomicroscopy was tedious, for in 
nearly all cases the structures consisted of a hard and brittle precipitate 
(B-manganese, for example) in a soft matrix of silver. Despite extreme 
care, a certain amount of pitting is present in most microsections. Over 
fifty alloys were so examined, most of which contained less than 50 at. % of 
aluminium, but for the purposes of this discussion, remarks will be confined 
to the silver-rich corner of the phase diagram. 

X-ray investigations were used to complement the micrographic work, 
and in addition to the normal cameras, involved the use of a Unicam 8150 
High Temperature evacuable powder camera. Samples of filings were 
sealed in vacuo in thin-walled silica capillaries. 

Finally, the identification of the non-equilibrium precipitates in the dilute 
alloys around the composition Ag;MnAl could only be undertaken with 
extraction techniques familiar in electron microscopy, and using here a 
Metropolitan-Vickers EM3 machine, which has facilities for electron 
diffraction. 


§3. THe PuHase Diacrams 


The major difference between the copper-based and silver-based Heusler 
alloys lies in the behaviour of the 8 (body-centred cubic) fields of the two 
binary systems, copper—aluminium and silver—aluminium, occurring around 
25at.% aluminium. In the former the £ phase is in principle stable after 
quenching from elevated temperatures, although in practice it partially 
transforms to other complex products (Gawranek et al. 1936). Theaddition 
of manganese to this system apparently suppresses this martensitic trans- 
formation, and produces the highly ordered structure characteristic of the 
Heusler alloys. In the silver-aluminium system, however, the f phase is 
only stable at elevated temperatures and decomposes completely, again by 
a martensitic transformation, into a close-packed hexagonal structure (C). 
The equilibrium structure (1) below 450°c is isomorphous with B-mangan- 
ese. One of the first objects of this investigation was therefore to see if the B 
phase could be retained at low temperatures by the addition of manganese. 
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To reduce the large amount of labour normally involved in ternary 
diagram studies, the problem was reduced to a consideration of structures 
normally present under the following conditions : 


(i) The section of the equilibrium ternary diagram at 400°o. This 


temperature was chosen as close to but above the Curie points of the alloys 
concerned. 


(ii) The section of the equilibrium diagram at 650°o. This is within the 
B phase field of the binary system, but could only be examined using high- 
temperature X-ray studies. 

(i) The structures arising after quenching in water from 650°c. These 
non-equilibrium structures are related to the development of the 
ferromagnetic properties. 

Fig. 2 


Je) 
AZOOPRE O = SINGLE PHASE 
© = TWO PHASE. 
@ = THREE PHASE. 


Approximate equilibrium diagram for the silver-manganese-aluminium 
system at 400°c. 


A feature of the investigation was the rapidity at which equilibrium was 
reached, even at 400°c. Annealing for 24 hours at that temperature gave 
structures which differed little from those annealed for 7 days, although the 
latter period was adopted in all cases. At 650°C, equilibrium was attained 
in the course of a few hours, and the specimens could be successfully 
heat-treated within the high temperature x-ray camera itself. 
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The main intent in obtaining the ternary diagrams was to identify the 
phases appearing in the ferromagnetic alloys, and no high accuracy 18 
claimed for the diagrams. However, they appear essentially simple, no 
intermetallic compounds have been discovered, and it is felt that the 
number of alloys examined enables the trend of the boundaries to be given 
with fair certainty, at least in the silver-rich region. Thus, only about half 
the alloys examined appear marked on the diagrams, the remainder (of 
lower silver content) merely served to confirm the run of the boundaries. 


Fig. 3 


650°C 


O = SINGLE PHASE. 
@© = TWO PHASE 
@ = THREE PHASE. 


Ag NG : sO 


Approximate equilibrium diagram for the system at 650°c. 


The diagrams are the combined results of x-ray and micrographic examina- 
tion in conditions (i) and (iii) outlined above, while the section at 650°c is 
from high-temperature x-ray investigations alone, with correspondingly 
lowered accuracy. 

Figure 2 therefore shows the 400°c section of the silver-rich corner of the 
equilibrium diagram for specimens quenched from that temperature. It 
will be seen that the alloys around the composition Ag;MnAl lie within a 
broad two-phase field where the equilibrium structure is a—Ag+$—Mn. 
The absence of intermetallic compounds show a further difference between 
this system and the Cu-Mn—Al series (West and Thomas 1956). The pu 
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phase, although isomorphous with 8-Mn, has only a very limited extent, 
and shows no tendency, as might be expected, to form a single phase band 
with its isomorphous structure in the manganese-aluminium system. 
The direction of this field may fairly reliably be determined by the con- 
stitution of the alloy Ag,,Al,Mn, where the minor constituent (w) could be 
detected by x-ray examination because of its high scattering factor. 
The extent of the ¢ field is not well established, but again is not of importance 
in this discussion. 
Fig. 4 


QUENCHED FROM O = SINGLE PHASE. 


650°C. © = TWO PHASE. 


@ = THREE PHASE. 


Ag 
Approximate constitution diagram for the system, quenched from 650°c. 


Figure 3 shows the section at 650°c. As in the copper-manganese— 
aluminium system, the f field has extended to cover a large area of the dia- 
gram, and overlaps, though incompletely, the area of the strongly ferro- 
magnetic alloys. A limited number of investigations at even higher 
temperatures (up to 790°C) showed that the B phase could extend still 
further, close to but not embracing the Heusler composition Ag,MnAl. 
It is therefore apparent that the strongly ferromagnetic alloys contain a 
most only a small proportion of 8 at 650°c—those around the ‘Silmanal 
composition contain none. 

Finally, fig. 4 shows the constitution diagram of alloys quenched in water 
from 650°c. Inno case was it found that the B phase was stabilized by the 
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addition of manganese, in striking contrast to the Heusler alloys. A large 
single-phase loop of Cis therefore formed, which, dependent on composition 
will generally decompose into a~Ag + B-Mn on reheating to elevated tem- 
peratures. The other point of note here is that the increased solubility 
of manganese in silver gives rise to an enlarged w solid-solution field, and 
precipitation of B—Mn is therefore to be expected in the ‘ Silmanal’ region at 
higher temperatures. 
§ 4. MicROSTRUCTURES 


The optical micrographic work as mentioned in § 2 was an integral part 
of the phase diagram determinations of the previous section, and helped 
to identify the structures present, particularly at low precipitate concen- 
trations. Certain new features however are relevant to the formation of 
the magnetic alloys. 


4.1. The 400°C Section 


The three-phase field (w-Ag ++ f-Mn), shown in the 400°c section, is 
easily identified ; an acid ferric chloride etch proving suitable here. Figure 
5 (Pl. 86) shows the alloy Ag,A1,Mn, which is close to the centre of this field. 

The two-phase a—Ag + B—Mn field within which the magnetic alloys lie is, 
however, moreimportant. Alloys low in silver show a clumped precipitate 
of B—Mn in the « matrix. 

Figure 6 (Pl. 86) shows the alloy Ag.MnAl, with a structure characteristic 
of alloys containing 50 at. % or lesssilver. As the silver content is increased, 
however, the alloy Ag,MnAl shows (fig. 7, Pl. 87) that a finely divided 
Widmanstitten precipitate has developed, while the amount of the 
clumped precipitate is diminishing. The proportion of the Widmanstatten 
pattern increases in the alloy Ag,MnAl, (fig. 8, Pl. 87), while in ‘ Silmanal’ 
itself Ag;MnAl (fig. 9, Pl. 88), the only clumped precipitate is a small amount 
left in the grain boundaries themselves, and all the grains are covered with 
the Widmanstitten precipitate. Thereafter, the amount of precipitate 
decreases rapidly as the « boundary is approached, so that in Ag,MnAl only 
a very fine precipitate is left, just within the limit of resolution of the 
microscope. 

This Widmanstatten pattern is therefore a feature which, in its develop- 
ment, follows the magnetization contours of Potter’s alloys (fig. 1), and 
therefore suggests it is the remains of a magnetic phase. However, it 
cannot constitute a ferromagnetic precipitate, since these micrographs 
(figs. 5-9) are taken from specimens annealed at 400°c, where, as Potter 
(1931) showed, the magnetization is rapidly destroyed. The structure 
of the precipitate below the Curie point is therefore of prime importance, 
since the equilibrium structure (8-Mn) is always paramagnetic. 


4,2. Non-equilibrium Precipitates 
Four specimens of ‘Silmanal’ (Ag;MnAl) were then heat-treated at 
790°C in an evacuated silica tube and quenched in water. One specimen 
(magnetic) was examined in this state, two others were annealed for 7 days 
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at 240°C and 310°C respectively yielding other ferromagnetic alloys, and 
one for 7 days at 400°c. The latter was tested and found non-magnetic. 
The specimens were then polished and carefully etched in cyanide—persul - 
phate solution, and examined under high power to observe the development 
of the precipitate. In the case of the first two magnetic alloys, (i.e. unaged 
and aged at 240°c), the precipitate was barely resolved. Figure 11 (PI. 89) 
shows the specimen aged at 310°c—the large specks are undissolved 
B-Mn, but it is apparent from the etching attack that a fine background 
precipitate has formed which must be the ferromagnetic phase. In the 
specimen aged at 400°, however (fig. 12, Pl. 89), the typical Widmanstiitten 
precipitate has again appeared. The higher ageing temperature should 
yield a coarser precipitate, but at this stage the transformation to B-Mn 
has taken place. It is therefore concluded that the ferromagnetic constitu- 
ent is an intermediate phase, formed during the nucleation of B-Mn, which 
can grow into platelets crystallographically related to the matrix, and which 
transforms to the equilibrium structure at temperatures above the Curie 
point. A final solution to this problem could therefore only be reached by 
extraction of this constituent itself. 


§ 5. ELecTRon Microscopy 


Extraction of the ferromagnetic precipitate was first attempted by 
electrochemica] means. A suitable specimen of ‘Silmanal’ was made the 
anode of a bath of 10% citric acid in water; the cathode was a silver foil. 
The sludge formed was collected on a sintered glass filter, dried, and further 
separation carried out over an alternating current electromagnet. Despite 
many attempts, the separation thus achieved always yielded a sample whose 
x-ray patterns were that of silver, and the technique was abandoned. 

A successful extraction technique was eventually developed using the 
carbon replica technique familiar in electron microscopy. A film of carbon, 
some 50-1004 thick, is evaporated onto the polished and lightly etched 
surface of the specimen, which is then immersed in a 10°% cyanide—persul - 
phate bath. The etch diffuses through the carbon film and attacks the 
matrix, and, after sufficient time has elapsed, the carbon can be floated off 
on water, carrying with it not only a replica of the surface, but also the 
precipitate particles embedded in it. Since this etchant is normally 
slow, prolonged immersion was usually necessary ; periods up to one hour 
were common. In obstinate cases, the films were finally removed by a 
brief immersion in 10°% nitric acid in alcohol. 

By this means, replicas of the magnetic and non-magnetic samples 
mentioned in §4 were obtained, from which high-resolution electron 
micrographs and also electron diffraction patterns could be obtained. 

Figure 13 (Pl. 90) shows the structure of the sample aged at 240°C. The 
particles are seen to be irregular in shape, and generally less than 10004 in 
diameter. A few examples of needle-like particles were seen, roughly 
50004 long and 500A across, which presumably were to become the 
Widmanstitten pattern appearing at higher ageing temperatures. The 
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electron diffraction pattern of the precipitate is shown in fig. 14 (Pl. 90). 
Not all the diffraction rings can be reproduced because of the high contrast, 
but on the original nine rings can be distinguished, and their Bragg spacings 
(d) determined. A small correction to the d values is needed, since the 
accelerating voltage is not exactly the nomina] 75kv. This can be made by 
evaporating a layer of thallium chloride (whose structure is known) onto the 
surface of the replica, and plotting a graph of true d values against measured 
d values. A straight line relation is obtained, whereby the d values of the 
precipitate can be similarly corrected. 


Analysis of Electron Diffraction Pattern from ‘Silmanal’ Extract 


Ring | Radius | Corrected d kil} 


no. (em) (A). 
i 0:53 4-72 {0001} 
2 | 0-94 2-66 1010} 
o. {0002} 
3 1-07 2-34 {1011} 
4 1-4] 1-78 {1012} 
{0003} 
5 1-59 1-58 11303 
ies {1013} 
6 | 1-84 1-36 \ pare 
.9§ ; {1122} 
a 1-92 1:31 20313 
aoe {0004} 
2 a {2022} 
9 2-33 1-07 {1014} 


The table shows the results for the extract. The analysis of the results 
is achieved by a consideration of the possible structures which may arise. 
The Heusler structure has certain characteristic absences, and the intense 
lines have hkl values such that }/?=8, 16, 24, 32... and weaker super- 
lattice lines with }A?=3, 4, 11, 12, 19.... Although the intensities of 
x-ray and electron diffraction patterns are not strictly comparable, the 
intense lines certainly do not follow this pattern, nor in fact can a good 
cubic parameter be found for the lines. On the other hand, it is easy to see 
that the d values of lines 1, 3, 5, and 8 bear a simple ratio to one other, 
and trial plotting with Bunn charts allowed all the lines to be satisfactorily 
indexed, giving a primitive hexagonal cell with parameters a=2-944, 
c/a=1-61. This phase will be designated ¢’, and must therefore be the 
ferromagnetic constituent. 

Replicas from the specimen aged at 400°C revealed the typical 
Widmanstatten structure which had been observed in figs. 9 and 12, but in 
greater detail (fig. 15, Pl. 91). Extraction of the coarser precipitate 
proved difficult here, but satisfactory evidence was obtained that the 
structure present in this state was predominantly B—-manganese (a = 6-39 A). 
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§ 6. Discussion 


The results so far described have shown that the presence of ferro- 
magnetism in the silver-based Heusler alloys is due to the presence of a 
non-equilibrium precipitate which forms on ageing from the super-saturated 
silver matrix, grows into a Widmanstiitten structure, and eventually 
transforms above the Curie point into the equilibrium product, p- 
manganese. 

Itis only possible at this stage to say little more about the crystal structure 
of this ferromagnetic phase. However, the existence of strong diffracted 
intensities of the form {000/} with / odd show that the lattice in its magnetic 
state is ordered. It is not, however, possible to distinguish between say, a 
nickel arsenide structure and any other hexagonal type which gives rise to 
ferromagnetism. The nickel arsenide structure can probably be rejected, 
since this cell is somewhat smaller than the equivalent structures MnAs 
and MnSb, bringing the manganese atoms too close together. If, on the 
other hand, the ordered structure results from sheets of manganese and 
aluminium atoms, separated by other sheets of silver atoms, the closest 
distance of approach of the manganese atoms could be equal to the c axis, 
4-704. This compares with a manganese separation of 4-18 4 in Cu,MnAl, 
4-37 in Cu,MnSn, and 4-39 in Cu,MnIn. 


6.1. Comparison with Previous Work 


It was mentioned in § 1 that the only previous work on the structure of 
these alloys was that of Geisler (1951), and it is now possible to make a 
detailed comparison with the present results. 

In the course of a study on the structures of permanent magnet alloys, 
Geisler examined microsections of Ag,MnAl, and noted the presence of a 
precipitate with a Widmanstatten structure. After a consideration of 
alloys around this composition, he concluded that ‘Silmanal’ was a simple 
precipitation hardening alloy, and that the precipitate itself had the cubic 
Heusler structure. The present work shows that this is not the case, and 
the presence of a Widmanstatten structure is not necessarily evidence of a 
ferromagnetic state. 

The other point of note concerns the orientation relation between the 
Widmanstatten structure and the matrix. From the microstructures, 
Geisler concluded that the platelets were aligned with the {001} planes of the 
matrix. The lattice parameters found here suggest an alternative orienta - 
tion. Mehl and Barrett (1931) noticed the appearance of just such a 
Widmanstitten pattern in a binary alloy of silver and aluminium containing 
20-2% silver, where the equilibrium structure is a precipitate of Cin a face- 
centred cubic matrix. X-ray examination showed that the orientation of 
the precipitate was 

(0001), //(111)m» 


[1120],,//[110],, - 
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In the present case, the lattice parameter of the terminal silver solid solu- 
tion (4-064 4) is very close to that of the aluminium solid solution studied in 
the binary alloy; the hexagonal parameters likewise correspond. Calcu- 
lations then show that if this orientation relationship were maintained here, 
fitting would be achieved with a linear strain of 2-3% along the a-axis, and 
only 0:6°% along the c-axis. 

It is tempting to pursue the analogy between the present case and the 
binary silver-aluminium alloy. Walker and Guinier (1953) on the basis of 
low-angle x-ray scattering experiments showed that at early stages of the 
precipitation of ¢ in the binary alloy, the particles were spherical, and it was 
only at the later stages of ageing that the particles assumed a platelet form. 
A consideration of figs. 13, 11 and 15 suggest that this may also hold in the 
ternary system. 

It must not be assumed, however, that there is any close relationship 
between the ¢ phase in the quenched alloys or in the equilibrium structure, 
and the ¢’ phase responsible for the ferro-magnetism. For example, the ¢ 
phase never forms a Widmanstiatten pattern in quenched alloys, and if it 
were responsible for the ferromagnetism, the maximum intensities of 
magnetisation would lie much closer to the composition Ag,Al than they 
in fact do. Further evidence against this association will be given below, 
and since all the magnetic alloys lie in the binary field a—Ag+f/—Mn, @ 
would appear to be simply a non-equilibrium nucleation phase in the 
precipitation of 8B-manganese from the supersaturated matrix. 


6.2. Magnetic Measurements 


A large amount of unpublished data has been obtained in this Department 
on the magnetic properties of ‘Silmanal’ and associated alloys, but in this 
discussion attention will be concentrated on two of interest, (i) the variation 
of magnetization with temperature, and (ii) the origin of the very high 
coercivities in these alloys. 

The variation of saturation magnetization with temperature provides 
further evidence for the distinction between the ¢’ (magnetic) and the 
¢ (non-magnetic) phases. Figure 16 shows the variation of intensity of 
magnetization (/) against temperature, for the alloy Ag,Al,Mn, which is the 
alloy closest to Ag;MnAl where the precipitate € can be detected in the 
quenched state. The curves show that in the chill-cast or quenched 
condition, J is low but after a heating and cooling cycle to a temperature 
above the Curie point, / is markedly increased. The interpretation here is 
that the a+ ¢ structure in the as-quenched state transforms to supersatur- 
ated a or even «+ ¢’ in the transition to the equilibrium « + 8—Mn structure, 
and hence on cooling a large number of the ferromagnetic nuclei of 
¢’ remain. 

In *Silmanal’ itself (fig. 17) the result is very different. Here the chill- 
cast or quenched specimens show the maximum value of J, and any pro- 
longed heat-treatment above the Curie point reduces this by a large factor, 
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the ¢’ phase transforming rapidly to the equilibrium B-Manganese. Since 
both these alloys lie in the binary equilibrium field of a-Ag+f8-—Mn, and 
since I for Ag,Al,Mn is less than J for Ag;MnAl, there must therefore be 
important differences in structure or ordering between the ¢ and ¢’ phases. 


Fig. 16 
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Variation of intensity of magnetization with temperature for Ag,Al,Mn. 


One can further predict that all alloys in the a—Ag + —Mn binary region 
will be similarly ferromagnetic. Even those low in silver will dissolve a 
small amount of B—Mn on heat-treatment. Some of this will precipitate as 
¢’ on quenching, and ageing will render the alloy ferromagnetic. As an 
example, the alloy AgAIMn, was found to be weakly ferromagnetic after 
suitable heat-treatment. 
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6.3. Origin of the Coercive Force 


Kittel (1949) has reviewed the origins of high coercivities for small, single 
domain particles, which are generally stated to be due to either (i) shape, (ii) 
magnetocrystalline anisotropy, or (iii) strain. In the silver-based Heusler 
alloys a finely dispersed magnetic phase is present, with particles under 
10004 which are in all probability of single domain size. Further, it has 
been shown that the coherency strains are so low that this effect may be 
ignored. 


Fig. 17 
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Variation of intensity of magnetization with temperature for Ag,MnAl. 
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In high coercivities due to shape anisotropy, single domain particles can 
produce a maximum theoretical coercivity H, given by 


H,= 2, 


where J, is the saturation magnetization for material in the form of long 
cylinders. The electron micrographs, however (fig. 13) show that the 
precipitate particles in the magnetic state are irregular in shape and do not 
possess these high axial ratios. It is only in specimens heat-treated at 
higher temperatures (fig. 15) which these axial ratios of 10 or more are found 
for the precipitate particles, and by this stage the specimen is non-magnetic. 
Furthermore, approaching this state the particles are almost certainly larger 
than single domain size, and as Meiklejohn (1953) has shown, the presence 
of only a small proportion of multidomain particles will appreciably lower 
the coercivity of a specimen. Alternatively, one may exclude shape aniso- 
tropy as a possible mechanism from an estimate of J, for the precipitate. 
I, for the alloy Ag;MnAl reaches a maximum value of 55-65, and if it is 
assumed that the particles occupy less than 5°% of the total volume (cf. 
fig. 12), 1, for the particles may have a maximum value of 1300 and 
H, > 8000, while values of H,~ 8300 0e have been reported in ‘Silmanal ’ 
(Hollis, private communication). 

By the process of elimination, it is therefore more likely that magneto- 
crystalline anisotropy is the predominant factor, giving a theoretical 
coercivity 

H,=0-64K/I, 


for particles distributed at random with respect to the field, where K is the 
anisotropy constant. Lacking knowledge of the value of K, we can only 
assume that its magnitude, in common with other uniaxial crystals, is high. 
More significantly, however, a plot of H, against temperature should fall off 
more rapidly than a plot of J, against temperature, since anisotropy is more 
readily affected by increased temperatures than the saturation magnetiza- 
tion. Thisis, in fact, found to be the case (Hollis, private communication). 


6.4. Tart VALENCY OF MANGANESE 


The f phase fields of the Cu—Al and Ag—Al alloys are both examples of 
intermetallic compounds having a 3:2 electron—atom ratio. Coles et al. 
(1949) have suggested that the Heusler structure (Cu,MnAl) itself belongs 
to this category, and for this to be so, the manganese atoms must be 
exerting a valency of unity. On the other hand, Hume-Rothery (1948) has 
shown, from a consideration of the form and direction of the B phase boun- 
dary in Cu-Mn-—Al alloys, that the manganese atoms would be divalent, 
and it would not therefore seem to be legitimate to consider the Heusler 
phase as belonging to this group of structurally analogous phases. A 
similar conclusion arises from this study of the Ag—Mn—Al alloys. The 
direction of the 8 phase boundary at 650°c (fig. 3) is again consistent with 
manganese acting asa divalent metal: furthermore, the alloy of composition 
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Ag,MnAl does not come within this 6 phase field, even very close to the 
solidus at 790°c, and so lends support to this view. 

The AgsAl (2) structure, however, shows a different aeen Although of 
only limited extent at 400°c (fig. 2), the tangential direction of the phase 
boundaries is consistent with the manganese exerting a valency of zero, 
and we thus have the example of manganese exerting two different valencies 
in structures which are both considered to be examples of a 3/2 electron 
compound. 
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ABSTRACT 


The absorption of 13-6 Mc/s sound waves has been measured in liquid 
helium I under pressures of up to 25 atm. The absorption varies considerably 
with the density of the liquid, but at all temperatures above about 3°K is 
adequately represented by the classical Stokes—Kirchhoff expression. 


§ 1. INTRODUCTION 


THE absorption of sound in a fluid is given by the classical expression 
due to Stokes (1845) and Kirchhoff (1868) namely 


Peay (A oe 
a (A+ 2n)+ 1)— a nl) 


where fis the frequency of the sound, v the velocity of sound, p the density 
of the liquid, 7 the coefficient of shear viscosity, A the coefficient of dilational 
viscosity (see for example Wilks (1958)), y the ratio of the specific heats, K 
the coefficient of thermal conductivity, and C,, the specific heat at constant 
pressure. The coefficient A is not readily identified with any property 
of the liquid which may be easily measured, but it follows from kinetic 
theory (Maxwell 1867) that for a perfect gas 


| eee Ve ake ee Be 1s SP 9 oo |) 


so that (1) becomes 
Qar*f 2 K 

c= = [7+ Lo sg gor, Mt. ase en) 
The thermal motion of most liquids is generally more complicated than. 
that of a perfect gas and this results in the absorption being considerably 
greater (see for example Markham ef al. (1951)). However, measurements 
on liquid helium I under the vapour pressure have shown that in the case 
of this simple liquid the absorption is given by expression (3) (Pellam and 
Squire 1947). The only other liquids known to behave in a similar way 
are liquid mercury (Ringo ef al. 1947), liquid argon, liquid nitrogen and 
liquid oxygen (Galt 1948). The present paper extends the earlier results 
on helium I by making measurements with the liquid under pressures of 
up to 25 atm, so that several of the parameters in eqn. (3) take up consider- 


ably different values. 


en eeetes ee 
+ Communicated by the Authors. 
+ Now at Department of Physics, Institute of Cancer Research, London. 


746 J. A. Newell and J. Wilks on the 


§ 2, THE MEASUREMENTS 


The apparatus and methods were similar to those used in experiments 
on liquid helium IT above 1°K, which we have described previously (Drans- 
feld et al. 1958). The measurements were made at a frequency of 13-6 Me/s, 
and the values of the absorption obtained are shown in the figure for various 
pressures over the liquid. After making due allowance for the difference 
in frequency these results are in satisfactory agreement with Pellam and 
Squire’s measurements at 15 Mc/s. The second column of table 1 gives 


0.5 SVP 

Ona 

0.3 

0.2 3.9 atm 

0.1 8.4 atm 
16-4 atm © 
24.7 atm 


2-0 3:0 4.0 
° 
Pas | 
The coefficient of absorption of 13-6 Mc/s sound waves in liquid helium I as a 
function of the pressure over the helium. 


Table 1. The Absorption of Sound in Liquid He I under the 
Saturated Vapour Pressure 


4-0 0-43 97 48 0-86 56 104 
3:5 0-34 110 49 0-65 50 Oe) 
3:0 0-25 102 51 0-38 4] 92 
2:5 0-20 87 48 0:07 6 54 


a is the measured absorption, «, is the reduced absorption, «,* is the reduced 
absorption calculated from (4) using values given in columns 4, 5 and 6. 
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smoothed values of the absorption « obtained under the vapour pressure, 
while the third column gives values of the reduced absorption «,, obtained 
by dividing the measured absorption by the factor 272 f?/pv®. Values of 
the density have been given by Keesom and Keesom (1933) and of the 
velocity of sound by Atkins and Stasior (1953). Values of « and x, at the 
higher pressures are given in the second and third columns of table 2. 


Table 2. The Absorption of Sound in Liquid He I under Pressure 


T°K - es 47 vy (y-1)K/C, ow, 
(em!) eos x 108 ” <0s x 108 
8-4 atm 
4-0 0-098 84 68 0:36 32 LOO 
3:5 0-087 85 70 0-26 26 96 
3-0 0-082 85 2 0-15 19 oI 
16-4 atm 
4-0 0-071 112 86 0-25 24 110 
3°5 0-064 106 89 0-15 ie 106 
3-0 0-061 102 92 0-10 14 106 
24-7 atm 
| 
4-0 0-050 119 102 0-18 18 120 
3:5 0-047 114 100 0-10 12 112 
3-0 0-046 ial iO 0-07 10 121 


aw is the measured absorption, «, is the reduced absorption, «,* is the reduced 
absorption calculated from (4) using values given in columns 4, 5 and 6. 


§ 3. DiscussIon 


It is convenient to re-write eqn. (3) using the reduced absorption «,, 
thus 


om Ite we oe Se Bee (4) 
To compare this expression with the measured absorption under the vapour 
pressure we use values of 7 given by Taylor and Dash (1957) and values 
of K given by Grenier (1951). C,, we take as being approximately equal to 
the specific heat under the saturated vapour pressure as measured by Hill 
and Lounasmaa (1957). Values of y and (y—1) are obtained from the 
compressibility measurements of Keesom and Keesom (1933) together 
with the velocity of sound. The fourth, fifth and sixth columns of table 1 
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give values of the various quantities on the right-hand side of eqn. (4), 
while the last column gives the calculated values of the reduced absorption 
(a,*). At 3°K and above the observed absorption is adequately accounted 
for by the classical expression, a conclusion previously reached by Pellam 
and Squire using earlier and less exact data for the various parameters. 
Any discrepancies are within the errors of the several experimental results 
used in the comparison. In particular, it is difficult to estimate values 
of y accurately for pressures close to the vapour pressure, also the velocity 
of sound comes into the expression to the third power. Below 3° the 
absorption is greater than is given by expression (4) and this is certainly 
due to an additional absorption associated with the Jambda transition 
which does not concern us here. (For further details see Pippard (1951) 
and Chase (1958).) 

The position at higher pressures is shown in table 2. The viscosity at 
higher densities has been measured by Tjerkstra (1952) and the third 
column of the table gives his values for 47. On comparing these figures 
with those for «,, itis seen that most of the absorption at the higher pressures 
is due to the viscosity. The experimental data available for evaluating 
the contribution to the absorption which arises from the thermal con- 
ductivity, is rather sparse, and we have to make the best estimate we 
can. The thermal conductivity under pressure has been measured by 
Grenier (1951) although only at 4:2°K and up to 3atm; there was no 
significant change with pressure. The absorption due to thermal con- 
ductivity is also proportional to (y—1)/C, or to (y—1)/yC,. The only 
measurements of the specific heat of liquid helium I under pressure are a 
set of values for C,, corresponding to a pressure of about 19 atm (Keesom 
and Keesom 1935), and these suggest that to a rough approximation C,, 
does not vary very much with pressure. Hence if we assume that AK and 
C, are sensibly constant the reduced absorption at any temperature will 
be proportional to (y—1)/y. Values of (y—1) taken from Atkins and 
Stasior (1953) are given in column 5 of the table, and estimated values of 
the thermal contribution to the absorption in column 6. Finally, the last 
column shows the total estimated reduced absorption («,**). Bearing in 
mind the uncertainties mentioned above which are involved in determining 
a, and «,*, there is reasonable accord between the estimated and observed 
values of the absorption. 


§ 4. CONCLUSION 


Above about 3-0°K the absorption of sound in liquid helium I at all 
pressures up to 25 atm is represented quite well by the Stokes—Kirchhoff 
expression (3). The large decrease in the absorption at the higher pressures 
(see figure) arises almost entirely from the large increase in the velocity 
of sound. Indeed the reduced absorption, «,, actually increases slightly. 
Although there is a rapid reduction in the reduced absorption associated. 
with the thermal conductivity, this is more than compensated for by the 
rise in the viscosity. 
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It appears that an essential condition for any liquid to have an absorption 
of the form given by (3) is that its molecules should have only translational 
modes of motion. Thus the liquids with such an absorption are either 
monatomic like helium, argon and mercury, or at so low a temperature 
that their rotational motion is frozen out as in nitrogen and oxygen. 
Although the Stokes relation, A= —3n, seems to be valid for these liquids 
it has only been proved for a perfect gas and there is no theoretical basis 
for applying it to liquids. Helium I under the vapour pressure has an 
extremely open structure, its viscosity measured at constant volume falls 
considerably with decreasing temperature (Tjerkstra 1952), and of all 
liquids approximates most closely to a gas. However, it is not clear that 
such an approximation is appropriate to other liquids, or even to helium I 
at the higher densities where the viscosity rises with falling temperature. 
A theoretical investigation of the value of the coefficient of dilational 
viscosity of simple liquids and dense gases would therefore be of interest. 
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ABSTRACT 


The results of palaeomagnetic measurements are given both for the 
Tertiary basalts of North Eastern Ireland, and for the Cretaceous sediments 
in the south of England. The direction of magnetization of the basalts is 
consistently shallower than the present dipole field direction by about 8°. 
The Cretaceous results are shown to be much less likely to be reliable, but 
also indicate a shallower dip than that of the present dipole field. A simple 
and reliable statistical treatment is developed for use with palaeomagnetic 
data. 


§ 1. INTRODUCTION 


Rock magnetic data for rocks in Britain are completely lacking for 
the interval between Tertiary and Triassic times, except for a few 
Jurassic measurements made by Nairn (1937). Since the Triassic results 
(Clegg et al. 1954) show a marked variation from the present dipole field 
direction while the later Tertiary rocks are magnetized, on the average, 
in a direction corresponding to that of an axial dipole, it is extremely 
desirable to gather information concerning intermediate times, in order 
to ascertain whether any progressive change took place during the interval. 
The author’s work has included two attempts to fill this gap in our know- 
ledge; one was a survey of the Antrim Basalts of Northern Ireland, and 
the other a series of measurements on Cretaceous sediments from the 
Isle of Wight. These investigations are described in this paper. 


§ 2. THe ANTRIM BasaLts 


During a survey of the basalts of North Eastern Ireland, the author 
has measured the palaeomagnetic directions of several igneous bodies. 
This part of his work constitutes an extension of Hospers and 
Charlesworth’s (1954) work in the same region. 

The igneous rocks in this region are of Tertiary age, generally accepted as 
Kocene, but there exists some evidence, based on a comparison of Scottish 
and Irish Fossil pollens (Hospers and Charlesworth 1954), that they may 
have been formed as late as Miocene times. They are basaltic, and may be 
subdivided into three categories for the purposes of this report; (A) the 
older ‘ lower basalts ’ (olivine), (B) the younger ‘ middle basalts’ (tholeiitic) 
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which were previously called ‘upper basalts’ (Patterson 1955), and (C) 
the intrusive bodies, dykes and sills, whose relative ages cannot always be 
determined. Grouping of palaeomagnetic measurements has been made 
on this basis. 

The principle underlying collection of samples was to collect as widely 
as possible geographically and temporally. The practical limitations of 
transportation often limited the collection to one sample per igneous body. 
Ninety samples were collected from fifty-seven, distinct igneous bodies, at 
forty-seven different sites extending around the north-east coast of 
Treland from Downhill, in the north-west corner of Co. Antrim, to the 
dykes exposed east of the Mourne Mountains. Collection was never 
made more than a few miles inland. 

The measuring procedure was to cast each sample of average mass 
6kg, oriented in a rectangular plaster block, and then to measure its 
direction of magnetization by means of the astatic magnetometer built 
by Blackett (1952). A core lin. in diameter and 13in. long was then 
cut from the block, and measured under a less sensitive magnetometer 
as a check on the first measurement. 

Table 1 shows the three mean directions derived from the writer’s 
measurements, and the equivalent data from Hospers and Charlesworth’s 
(1954) paper. JD is the magnetic declination measured in degrees east of 
North, J the inclination, . is the number of separate igneous bodies 
involved, NV is the total number of samples taken from the .” igneous 
bodies, and n is the total number of independent measurements made 
(plaster blocks and cores). 

It has been customary among rock magnetists to use a form of statistics 
developed by Fisher (1953). These statistics are based on an assumed 
Gaussian distribution of density of measured vector directions about the 
mean direction. In certain cases this distribution is found experimentally 
to be invalid, so that statistics based on any assumed density distribution 
must be held in suspicion until the distribution be examined experi- 
mentally. In particular, any ‘circle of confidence’ containing the mean 
direction with a specified probability may contain falsely implied precision. 
For this reason a simple system of statistical evaluation of rock magnetic 
data has been developed (Appendix) and used on the data in this paper. 
The three significant quantities used are the mean direction, the analogue 
on a spherical surface of the standard deviation for a single reading, 6, 
and the standard deviation of the mean, «. These can be found either for 
a set of igneous bodies, or for a set of measurements on a single igneous 
body. They are easily understood, and free from any assumptions about 
distribution. 

The measurements for a single igneous body were rejected whenever 
8 exceeded 10°. The directions were corrected for geological tilt, which 
in general was small. In cases where only one sample has been taken 
(N=1, n=2) from an igneous body, 6 has been redefined as the angle 
between, the two directions measured for that sample. 
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It can be seen from table 1 that the four results are reasonably consistent, 
all reversed, and that all of them indicate a deviation from the present 
dipole field of about 8-4° towards shallower inclinations. This is a signifi- 
cant amount in terms of e. 

It is necessary to consider the possibility that this deviation might be 
due to a magnetostrictive effect caused by the pressure of the overlying 
lavas. The geology of the region indicates that the lowest of the basalts 
have not been buried more than about 600 metres deep. The pressure 
at this depth is not great, and similar pressures applied in the laboratory 
made no discernibl. directional change. There is no obvious discrepancy 
between the inclinations of the middle and lower basalts, even though they 
must have undergone quite different pressures during the geological time 
involved. 


Table 1 
D(ETN) rf NE EN n R ) € 
| Lower basalts ilies: — 64-6° up 19 | 35 | FO | E7500) 22 hea 

Middle basalts 206° — 62-2° up 6 Ki 14 Dei 2a, Veli ee 
Intrusives 184° — 63-3° up 16 | 25 | 49] 15-268 | 17° | 4-4° 
Above three 

combined 183° — 64:1° up 41 | 67 | 133 | 38-45 ZO Se 
Hospers’ results 194° — 60-2° up 24 | 72 | 120 | 23-253 | 14°.) 2-0° 
Four results 

combined 1872 — 62-4° up 65 | 139 | 253 | 61-63 1s aie pc 
Present dipole 

field Oe +70-8° down} — | — | — — — | — 
Present field 347° +70-0° down | — | — 


(For definitions of symbols, see text.) 


On the assumption that the axis of the magnetic dipole field coincides 
on. the average with the earth’s rotational axis, the interpretation of this 
8-4° deviation as being due to land-mass movement would indicate that 
Ireland has moved northwards 11-5° of latitude or about 800 miles 
since the eruption of the basalts. 


§ 3. THe Ise or WIGHT CRETACEOUS SEDIMENTS 


There having been no igneous activity during Cretaceous times in 
Britain, a systematic attack was made on Cretaceous sediments from 
eight different localities in Southern England. Forty-six samples were 
collected from Compton Bay, Whitecliff Bay, Brixton Bay, and Sandown 
Bay on the Isle of Wight, and from Swanage, Durdle Door, Dungy Head, 
and Lulworth Cove along the Dorset coast. Each site presented cliff 
exposures. 

Apparent magnetic stability was found in only nineteen of these samples, 
all of which were collected from the four miles of continuous Wealden 
sediment exposure along Compton and Brixton Bays. These nineteen 
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samples consisted of buff sandstones and red and grey clays. ‘The other 
twenty-seven samples were rejected on the basis of the laboratory and 
geological tests mentioned below. From the nineteen samples, ninety - 
seven oriented 2cm cubes were cut and their directions of magnetization 
determined. The mean intensities were 1 x 10-6 and 5x 10-6 e.m.u. /em3 
for the sandstones and clays respectively. Table 2 gives the mean direction 
and standard deviation of the nineteen samples with respect both to the 
bedding planes and to the present horizontal plane. It also shows the 


Table 2 
D(ETN) | J(down) ny) € 
Present dipole field 0° + 67-6° — = 
Present earth’s field 351° + 64-0° — — 
Mean direction with respect to 
horizontal plane 345° + 63° 4-9° | 1-1° (N=19) 
Mean direction with respect to 
bedding planes Pas +59° 85° | 2:0° (N=19) 
Reverse of Eocene field direction 
(for Northern Ireland) re + 62° 18-5° | 2:3° (N =65) 
Triassic field direction (for Great 
Britain) 34° + 26° — — 


directions of the present earth’s field and the dipole field at the Isle 
of Wight, and the Eocene and Triassic field directions for Britain. 
Evidence concerning the magnetic stability of these nineteen samples 
was as follows: 


(a) The mean direction of magnetization with respect to the present 
horizontal plane lies only 4° (3-5 standard deviations of the mean) from 
the direction of the present earth’s field. 


(b) The slight folding of these Wealden beds allowed us to calculate 
that the 8-5° standard deviation with respect to the bedding planes reduced 
to 4:9° with respect to the horizontal plane. This conforms with Graham’s 
(1949) test of stability for folded strata. 


(c) Alternating field demagnetization of the sandstones, in zero steady 
field, revealed a soft magnetic component, which disappeared for peak 
alternating fields less than 180 oersteds and had an initial magnitude 
which was small compared with a harder component. This harder 
component remained constant in magnitude and direction until 180 peak 
oersteds was reached, where it began to be unpredictably affected. 
Removal of the soft component improved the agreement between measure- 
ments on cubes from any one sample. 


(2) Measurement of the magnetic moment of sandstone cubes suitably 
oriented in the earth’s field revealed a soft component, which grew 


P.M. 30 
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exponentially with time, and was small in magnitude with respect to a 
harder component. The harder component was stable over at least several 
months. 


(ce) Two superficially different types of rock, with different natural 
intensities, comprised the nineteen samples. 


The fact that the mean direction of magnetization is roughly 
intermediate between the Eocene and Triassic directions cannot be 
taken as proof of stability. Test (a) admits the possibility of the harder 
component possessing instability over a time greater than several months, 
but small with respect to the times involved in secular variations. It 
could, on the other hand, be a chance result of the particular tilt of these 
Wealden beds. Test (b) indicates instability of the harder component. 
Tests (c) and (d) reveal only that relatively soft and hard components do. 
exist, and tell us nothing of magnetic stability over geological times, 
except that the alternating field test is consistent with the interpretation 
of the harder component as being thermoremanence in magnetite, following 
the work of Thellier and Rimbert (1955). Fact (e) is circumstantial 
evidence in favour of stability. 

There is, therefore, doubt as to the validity of any palaeomagnetic 
interpretation of these results, and further work, preferably on similar 
sediments differently tilted, must be done to establish a conclusive result. 
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APPENDIX 


A statistical system has been used which employs the analogues, in 
the case of distribution of points on a spherical surface, of the standard. 
deviations used for a distribution of points on a plane. 
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(a) ‘The mean direction of N directions is found in the customary way. 
if A, Hy v be the direction cosines of the mean direction and A,, p;, v; the 
direction cosines of the NV measured directions, then A, 1, v are defined as 


A=SAIR, w=LwlR, v=DvilR 
where R= V/{(2A)? + (Sm)? + (D%)}- 
That is, the mean direction cosines define the mean direction. 


(6) The mean cosine of the deviations 8; from the mean direction is 
found in order to establish the standard deviation for a aie reading. 


PES ie > cos 8; = Benes ieee = te 
N 


whence the deviation d=cos!(R eh y 


This is the analogue of the ‘standard’ or r.m.s. deviation, since for 
small 6 
cos§ = 1—182=R/N. 
This gives the r.m.s. deviation 6, = «/{2(N — B)/N}. 


T.m.s. * 


(c) Then by analogy with the case of a planar distribution the standard 
deviation of the mean direction, «, may be taken as 


e=6//N. 
The expressions § and « are simple to use and understand, and express 
the facts adequately without recourse to assumptions about distribution 


of vector densities. They are also free of mathematical approximations 
which might limit their region of usefulness. 
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ABSTRACT 


An account is given of an experiment in which the thickness of the ° satu- 
rated’ helium film formed from the vapour on a stainless steel mirror was 
compared with the thickness of the film formed from the liquid phase. The 
thickness of the film formed from the vapour was almost the same above 
and below the A-point and of the same thickness as the helium I film formed 
when the mirror was in contact with the liquid. The thickness of the 
helium IT film, with the mirror in contact with the liquid, was much thicker 
because of the more stringent isothermal conditions prevailing on the mirror 
surface. 


§ 1. INTRODUCTION 


FRENKEL (1940) and Schiff (1941) have suggested that the formation 
of the helium film is essentially due to the van der Waals adsorption 
forces between the helium atoms and those of the substrate. Minimizing 
the sum of the van der Waals potential energy and the gravitational energy 
of an atom at the surface of the film leads to the following equation for the 
film profile : 
d=k/ Hs 

where £ is a constant and d the thickness of the film at a height H. This 
equation has been derived in another way by Rice and Widom (1953) 
and independently by Atkins (1954) and Meyer (1955). They consider 
the helium film as an adsorbed layer in equilibrium with a very slightly 
unsaturated vapour. An interesting feature of this idea is the suggestion 
that the profile and thickness of the ‘ saturated’ film should, for the strictly 
isothermal case, be the same whether the substrate is in contact with the 
bulk liquid or not, apart from effects due to curvature at the edges of the 
substrate (Meyer 1955). The purpose of the work described in this paper 
was to compare qualitatively the thickness of the film formed from the 
vapour with that formed when the mirror was in contact with the liquid 
surface. 

Due to the mobile helium I film it was essential, if the present experiment 
was to provide unambiguous results, to measure the thickness of a film 
formed on a surface which had no mechanical contact with the bath 
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liquid or the film arising from it. To achieve this requirement a mirror 
was suspended in the vapour phase from a small cylindrical magnet which 
floated above a superconducting lead dish. One advantage of such a 
system is that the heat input to the surface is limited to radiation only. 
The disadvantage is that, at temperatures below the A-point, the surface 
is only strictly ‘isothermal’ with the bath temperature when the mirror 
is in contact with the liquid bath. 

For an account of the various theories of the film and the results for 
the film thickness obtained by many authors, the reader is referred to the 
review article by Jackson and Grimes (1958). 


§ 2. THe APPARATUS AND MetTHop 


The apparatus used was essentially the same as that described in an 
earlier paper (Ham and Jackson 1957: see also Grimes and Jackson, to be 
published shortly in the Philosophical Magazine, referred to as paper II). 
The only essential difference between the present arrangement and that 
described in paper II is in the design of the experimental chamber (the 
mirror tube); this is shown in fig. 1. The lead dish C rested on the glass 
rim of the mirror tube and had a small hole through its centre. A copper 
tube T was attached to the lead dish and made central in the mirror tube 
by suitably spaced phosphor bronze strips. The upper part of the mirror 
tube terminated in a kovar-glass seal S and this was joined by Hoyt metal 
to a kovar cap K. This was joined to a length of german silver tubing D 
which suspended the mirror tube in the cryostat and also served as a 
connecting line to the reservoirs of helium gas which was condensed during 
an experiment into the mirror tube. The small permanent magnet M 
was supported on a brass stirrup which was connected to the top of the 
Dewar vessel by a length of german silver tubing coaxial with D. The 
Dewar head was designed to provide a means of rotating and lowering or 
raising the magnet to any desired position in a vacuum of 10~7 mm Hg. 
The stainless steel mirror G was 1 cm long, 4mm wide and 2 mm thick 
and was suspended from the magnet by a thin tungsten rod R. 

The pre-run treatment and method of collecting helium gas from the 
surface of the liquid in the cryostat and condensing it into the mirror 
tube was similar to that described in Ham and Jackson (1957) and paper IT. 
The principles of the optical method used for the measurements have been 
described in previous papers (Burge and Jackson 1951, Ham and Jackson 
1957) and will not be described here. 

Because of the poor thermal contact between the contents of the mirror 
tube and the liquid helium bath it was necessary to have a slight pressure 
of helium gas in the mirror tube during the cooling period, prior to trans- 
ferring liquid helium into the cryostat. The resulting pressure at the 
liquid helium temperature was always between 30 and 60% of the saturated 
vapour pressure. The nicol positions for intensity matching across the 
1/3 boundary of the barium stearate layers on the mirror were taken under 
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these conditions and used as the zero position of the apparatus; due 
allowance was made for the birefringence effects in the glass of the mirror 
tube caused by the pressure difference across the glass surfaces. This 
procedure was considered justifiable as an experiment of Bowers (1953) 
showed that about two or three helium layers are adsorbed under the above 
saturations and this is equivalent to a nicol change of about 3 min of arc. 


Fig. 1 


The mirror tube. 


The magnet, during the cooling of the apparatus, was suspended by 
the brass stirrup well above the lead dish. Following the transfer of 
liquid helium into the apparatus the stirrup was lowered slowly until the 
magnet released itself and then lowered to the bottom of the lead dish. 
Floating the magnet in the above manner always gave stable floating 
positions and the frozen-in currents in the superconducting lead dish 
probably helped to stabilize the magnet. The floating position of the 
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magnet was not exactly reproducible and the zero position varied from 
experiment to experiment; this was not a disadvantage, however, for the 
present purpose. 

The temperature was determined by measuring the vapour pressure of 
the liquid on either an oil or mercury manometer and using the T 


L55 
vapour pressure tables (van Dijk and Durieux 1958). 


§ 3. THe EXPERIMENTS WITH THE FLoatiIna MAGNET 


At a temperature of 2-05°K, } litre of helium gas was condensed slowly 
into the mirror tube. The magnet remained stable and the initial nicol 
measurements indicated that a thick film was present on the mirror. 
It was first thought that the helium gas supply contained an impurity 
sufficient to affect the film thickness but later measurements showed that 
the initial matching position of 68° 36’ had changed to 67° 16’ and remained 
constant at this value. Subsequent condensations indicated that the 
initial high reading of 68° 36’ was due to a drop of liquid helium remaining 
on the bottom of the mirror. During later condensations the drop 
remaining on the mirror surface was evaporated by switching on the room 
lights: thickness measurements (nicol readings) were made after sufficient 
time had elapsed for the attainment of temperature equilibrium. After 
each condensation the position of the liquid level relative to the bottom 
of the mirror was recorded. The results indicated that the thickness 
of the film formed from the vapour phase remained practically independent 
of the position of the liquid surface. Condensations were continued 
until the liquid reached the mirror surface and immediately large nicol 
readings were obtained which indicated the presence of a much thicker 
film on the mirror surface. At a height of 0-06 cm the nicol reading was 
69° 0’. The temperature was then allowed to rise through the A-point 
to a temperature of 2-48°K. When temperature equilibrium had been 
attained the nicol readings were recorded and the mean value was similar 
to those obtained when the mirror was in the vapour phase at a temperature 
below the A-point. A second experiment was carried out at 2-03°K and 
similar results were obtained; in the vapour phase the thickness was 
practically independent of the position of the liquid helium level until 
the latter reached the mirror surface when the thick helium II film formed. 
After the completion of this experiment the temperature was raised to 
2-45°K and the mean nicol reading was again similar to that obtained 
when the mirror was in the vapour phase at 2-03°K. 

Two experiments were carried out at a temperature of 2-47°K. The 
times required for the attainment of temperature equilibrium were much 
longer but the experimental procedure was as described above. The results 
were very similar to those obtained at temperatures below the A-point. 
The only important difference was that the thickness of the film formed 
on the mirror, when the liquid was in contact with the mirror was practi- 
cally the same as that obtained when the mirror was in the vapour phase. 
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As the magnet floated in slightly different positions during each experi- 
ment the results have been normalized to the results at 2-47°K, so that 
all the nicol positions agree at a height of 3cm. These normalized results 
are plotted in fig. 2. The sensitivity of the apparatus was not determined 
directly but estimated using the thickness measurements of Ham and 
Jackson (1957) for the height range 0 to 2:5 cm, at temperatures below 
the A-point. The estimates of the film thickness are therefore somewhat 
approximate, but the results show clearly that a film forms the vapour 
phase at temperatures above and below the A-point and has a thickness 
similar to a helium I film formed on a surface which is in direct contact 
with the liquid helium. 


Fig. 2 
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The thickness of the helium film formed (a) from the vapour phase (above 
0-3.cm) ; (b) from the liquid phase (below 0:3 cm). 


§ 4. Discusston 


The apparent discontinuity in the curve of nicol positions against 
height when the level of helium liquid reaches the mirror, obtained at 
temperatures below the A-point, seems to lend support to the idea that the 
presence of the helium film has something to do with the A-point pheno- 
menon of the bulk liquid. This would only be true if the temperature of 
the mirror surface, when in the vapour phase, was the same as that of the 
bath liquid. However, observations showed that the mirror was at a 
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shghtly higher temperature than the bulk liquid and this alters the inter- 
pretation of the results. Part of the energy of the light beam is absorbed 
at the surface of the mirror. 

Further evidence of this heat input was obtained at temperatures above 
the A-point when the mirror was in direct contact with the bulk liquid He I. 
A fine stream of bubbles of helium gas was observed to rise from a small 
hole in the bottom of the mirror. From an estimate of the size of the 
bubbles and their rate of formation the heat flowing to the liquid through 
the bottom of the mirror was estimated to be of the order of 10-° cals/sec. 
As the surface area of the mirror was of the order of 10 times that of the 
bottom of the mirror a nearer estimate of the total heat input would 
probably be 10 eals/sec. 


Fig. 3 
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When the mirror is surrounded by the vapour phase the heat absorbed 
at the surface is removed through the vapour phase (see paper IT): the 
average energy of an atom leaving the film is 3/2 kAT greater than a 
helium atom entering the film. This temperature difference A7' is present 
at all temperatures when the mirror is in the vapour phase and is also present 
when the mirror is in contact with the liquid helium at temperatures above 
the A-point. At temperatures below the A-point, when the mirror makes 
contact with the bulk liquid, the superfluid destroys any AT present 
and the film is effectively at the same temperature as the liquid helium 
bath. The presence of a AT’ on the mirror produces drastic effects on the 
helium film profile and has been described fully by Ham and J ackson 
(1957) and Jackson and Grimes (1958). Under this condition the thickness 
measured at a height H becomes equal to the thickness of the ideal 
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isothermal film at a height (H + AH) where AZ is proportional to AT’. ‘The 
effect of AT on the film thickness is shown in fig. 3 at a temperature of 2-5°K. 
The values of AH which fit the estimated results for the thickness of the 
helium I film are in the range 100 to 200 and the corresponding values of AT’ 
at 2:5°K are between 10-8 and 2 x 10-3 degk. The value of A7’ calculated 
from the results of paper II was 10-4 deg xk. As the mass and surface area 
of the present mirror are very much smaller than that of the mirror used in 
paper II the temperature excess is expected to be greater for the present 
mirror. 

As the AT exists at the helium II temperatures when the mirror is 
in the vapour phase we should expect the thickness of the helium II film to 
be approximately equal to the helium I film: this is in agreement with the 
results. This would not be the case when the superfluid helium came in 
contact with the mirror. One would expect the temperature difference 
to be somewhat greater when the mirror is in the vapour phase than when 
it is dipping into the liquid He I and anticipate, as a result, a small dis- 
continuity in the helium I results when the liquid level reaches the mirror. 
The scatter of the results prevented the resolution of this discontinuity 
if present. 

Finally it may be concluded that if the thickness of the helium I film 
could be measured under isothermal conditions, the resulting thickness 
would be of the order of that of the helium II film. Whether it would be 
thicker than the helium II film just above the A-point, as Franchetti 
(1957) suggests cannot be established from the present observations. The 
results are in qualitative agreement with the earlier results of Ham and 
Jackson (1957) and Grimes and Jackson (to be published) and indicate 
the correctness of the theories of the helium film which are based primarily 
on the van der Waals attraction. 
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ABSTRACT 


A sealar wave with random variations of amplitude and phase across the 
wave-front is assumed as a simple model of a radio wave after it has left the 
ionosphere. The second-order moments of the in-phase and quadrature 
components of the fluctuating part of the field are found using the Fresnel 
diffraction formula. It is assumed that these two components have a joint 
Gaussian distribution, and a parameter py, called the intrinsic correlation 
of the fading, is introduced ; this serves as a measure of the eccentricity of the 
ellipse of constant probability. It is shown that py is the amplitude of the 
diffraction pattern produced by a distribution of intensity proportional 
to the spatial correlation function of the fluctuations over the wave-front. 
As the wave propagates the correlation tends to zero. It is shown that for a 
given irregularity size and wavelength, the irregularities may be supposed 
to be situated at any height, up to a certain maximum height, at which 

~*" py=l. An analysis of day time fading on 16ke/s has been made. It 
appears that the intrinsic correlation of the wave at the ionosphere is very 
near one, and that the ionosphere imposes phase modulation on the wave. 


1. INTRODUCTION 


THE application of the theories of statistics and random noise in radio 
propagation has progressed steadily since the paper by Ratcliffe (1948) 
where it was first suggested that a fading radio signal was essentially the 
same as electrical noise that has been passed through a narrow-band filter. 
The two well-known papers by Rice (1944) are devoted to the study of 
electrical noise in circuits. In such theories of electrical noise and fading, 
the random part of the radio frequency current is assumed to be randomly 
phased (i.e. that the phase is uniformly distributed from 0 to 27), and it is 
necessary to consider only its average power spectrum. Now when one 
considers thermal noise or thermal radiation, which originates in the 
independent motion of a collection of electrons, it is clear that the total 
current or electric field is randomly phased. However, the fading of a 
radio wave is quite different in this respect ; for the initial source is the 
coherent radiation of a transmitter, certain modifications in amplitude and 
phase being introduced by the ionosphere when it is not perfectly specular. 
It is not obvious that the wave is randomly phased in this case. 

It is the purpose of the present paper to examine in detail how the wave 
becomes randomly phased as it propagates away from the ionosphere. 


eer oe ee 
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It will be shown how one may infer from the statistics of amplitude and 
phase at the receiver the sort of modulation that was imposed on the 
reflected wave and the size of the irregularities in the ionosphere. It is. 
usual to measure the irregularity size by correlating the fading at two 
spaced receivers. In the present method the signals which are correlated 
are from the same receiver but have a 7/2 difference of phase. 

In two recent publications Bowhill (1956, 1957) calculated the statistical 
properties of a fading signal. His results are similar to those given in 
$2 of this paper, but his derivation was based on the theory of the angular 
spectrum. The present method using Fresnel’s diffraction formula is 
rather simpler. 

Further information about the diffracted signal is obtained by assuming 
that the in-phase and quadrature-components of the signal have a joint 
Gaussian distribution. Instead of Bowhill’s correlation of amplitude and 
phase we shall define a parameter called the intrinsic correlation of the fading 
signal. It will be shown that this parameter is related to the diffraction 
pattern due to a single irregularity, so that it is then possible to give a simple 
physical interpretation of the changes in the signal as it propagates. 


2. DIFFRACTION THEORY 


We have to consider the propagation of a plane wave after it has been 
modified by irregularities which alter its phase and amplitude. In the 
following work the Fresnel approximation is used in calculating the signal 
at the receiver. ‘The complex amplitude is assumed known over an infinite 
plane, the direction of propagation being normal to this plane. The zy 
plane, which is called the screen, contains the known amplitude field. 
The receiver is at (21, y,, 2) and the propagation is in the +z direction. 
We shall later substitute a wave diverging from a point source for the 
plane wave. 

Fresnel’s diffraction formula (Baker and Copson 1950), valid for r>A, is 


x, y= Qa 
E(z,,Y;; 2 mal]. »Y¥)exp ( tat — son) eee Ea dy (1) 
: is * A Ar 
where r= 4/ [22 yf, y— %) ys=angle between the direction of 


propagation, oe fe te joining oy y) and (x,, y,,2), A=wavelength, 
w=angular frequency = 27c/A, F =the complex amplitude of the wave at 
the screen, = F(z, y)+1F (a, y). 

We shall assume that the variation in the angle of arrival of the diffracted 
radiation is small enough that one may use the expansion of r in terms of 
(v—2x,), (y—y,)?, which is an approximation used by Fresnel. Since 
the width of the angular spectrum is determined by the gradient of phase 
of F over the screen, there is a certain restriction on this. For example, if 
the period of phase fluctuation over the screen is one wavelength, the 
amplitude of the periodic variation should not exceed one radian. It 
may in fact be shown that any field in the screen of shorter spatial period 
gives field rise only to evanescent waves which disappear within a few 
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wavelengths of the screen (Booker and Clemmow 1950). We have 
therefore 


E(x, Yy, 2 a: = { ies F(x, y) exp (=r [(w—a,)?+ (y—n)"1) dx dy 
x exp iwt— =) 
A 
on 
= [H,(x,, y,, 2) +tH,(z,, y;, 2)] exp (ior =). aa KE# 


The quantities F,, H, and F,, H, are respectively the amplitudes 
of the in-phase and quadrature components of the signal. These are 
equivalent to the amplitude and phase, but are much easier .to discuss in 
the analysis. We wish to make applications to experimental work in which 
both phase and amplitude are recorded, and not amplitude alone. 

This is the same diffraction formula that has been used by Booker et al. 
(1950), Ratcliffe (1948), and Bowhill (1956). They derived it from an 
exact diffraction theory based on the concept of the angular spectrum of 
the diffracted radiation (Booker and Clemmow 1950). Hewish (1951) has 
made the exact calculation in a case of a periodic screen where the period 
of the phase variation was 50 A with an amplitude of up to 4 radians; the 
angular spectrum was narrow enough that the Fresnel approximation 
might have been used. 

In order to study the statistical properties of H and F when the 
ionosphere is irregular, we consider a hypothetical assembly of screens— 
the parent population in the usual terminology of statistics. Averages 
will be with respect to this population. In practice one samples by taking 
arecord of fading at the receiver for as long as the fading remains statistically 
stationary in time, so that the sample is homogeneous. We denote by a 
bar over a quantity its assembly average. We shall use (2) to find the 
following quantities 


Wee, Les Wey et eee ee oe eee) 
We assume that F(x, y) is a stationary function over the entire zy plane. 
The probability distributions of H and F are discussed later. 
The following results are taken without proof 


comoo) 2 2 
{] exp (ilalat + y2) + be + ey + d)) da-dy = “exp (i | a- = ]): 
Be ee (4) 
were [Je & (+0) Audi =0(E 7). pees 


where a, 6, c, d are real, and 4(, 7) is the Dirac delta function. (4) may be 
derived from the Fresnel integral 


|. exp (iz?) dxe=4(1+1)1/7/2. 


) 
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The average of EL, H, is obtained from (2), dropping the 


Qriz 
exp: { = aoe +tw ‘ 


B= = {[- exp (— 2 [(e—a)*+ m1) de dy= F (6) 


—-O— 0 
Thus the average signal remains unchanged in course of propagation. 


We shall now find #2 and HE* from which one finds H,?, ete. The 
square of the right-hand side of (2) is expressed as a double integral, 


he Ba — =| | | i F(a, y) F(x, y’) exp (= 


x {(%—%,)? +G—wP te — @)+ (y'—y)?f ) da dy da’ dy’. 
Write a’ =a2+& y'=y+yn. Since F is a stationary random function 


F(a, y)Fl@+é y+ d , a pease of €,7 only. The wand y integration may 
be performed rom 


ee 2 
He w= 5 Tie NF e+ Bara exp | EE eee 


E? is therefore independent of (x,, y,) and His a stationary ee Jace 
EE* is also expressed as a double integral : 


0 60 0D 0D Fa mt 
Wi-s( {Fewer Fy) exp (F 


X {= (2)? +(e a)" —(y-y,)? + (y’ — yy) dx dy dx’ dy’, 


mu = ( {| [ P(x, y)F* (a+, y+n) ex(F 


x {— 2vé — 2yn ») dx a dé dy 


= hee F(a, y)F*(a7+&, y+n)a(é, n) dé dy 


BR FR. ie. H2+H2=F 2+ FY. eh: 
The average power is conserved during the propagation. From (7), (8) we 
may find the averages H,?, H,?, HE, : 
By? = 4(E + H*)? = 1(H + BP + 2BB*), 
E> -1(E—P* P= — + Pome), 
1 


BE, = — (E+ EE pa Ale — H*), 
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A is assumed that p,, py, py, 0 as £, n> ©. Clearly p}=p,=p;.=1 at 
= n — O. 

For brevity the integrals occurring in (7) when the spatial correlation 
functions are substituted, are denoted as follows: 


yd Poste - 
B,+tA\= se}. _mléa) exp (= (+9?) ) dé dn, | 
s il 00 60 —71 
B,+iA,= | i p2(é, 0) exp (= (S+9?)) dé dy, (10) 


. 1 
By +tAy= De | 


The final result giving the averages and variances of H,, E, are obtained 
from (7)-(10): 


_ Pao) exp (GE (E+ nt) db dy 


ae 2 


E=F,, H,=F;, 
Hy? — Hy =3(F?— F,?)(1—A,)+3(F,?—F2)(1+ A.) — (FF, —F Py) By», 
oH P= 4(Fy?— Fy2)(1+.A,) + 3(F2— F,2)\(1—A) + (FP, FP) By, 
LE, —H,E,= $(F 72-2 *)B, = +(F,?— 2”) By — (FP, PLP )Ayp. 
(11) 
As an example let the correlation functions be identical and of the form 
i 
ri=Ps=Pu=exp| — e+) | ee eee i 2) 
0 
Then we have 
py. eee Z 
NS a rery Leo = Tee 


where 
= arg | 2Az. 


This example will be discussed in § 3. 


2.1. Point Source 


As nothing has been said about what happens to the wave in the 
ionosphere, we cannot state generally what happens when the source is 
brought to a finite distance. Despite conflicting theories of the details 
of the scattering it has been customary to assume that the function F(x, y) 
is modified as follows for a point source: 


F : Q70 
Pex, y= (P(e y+ iF fe, 2} exp (iat a 


where 7, is the distance from some source on the —z axis at (z,, 0) to the 
point (x, y). a@=some constant. 

This is true only when the scattering occurs mainly in a region where the 
refractive index pis ~ 1. If there is appreciable scattering near the level 
of reflection—analogous to deviative absorption—then the modulation / 
just below the ionosphere will depend on, the detailed boundary conditions 
within the ionosphere (height gradient of electron density, etc.). It is a 


768 — R. P. Mercier on the 


mistake to assume always, for example, that the diffraction pattern on the 
ground is doubled in size when a point source is substituted for an incident 
plane wave. What follows therefore is correct only for a model of the 
ionosphere in which the scattering occurs in a region where p+ 1. 

We shall neglect the variation in intensity of illumination over the 
screen, and consider only the phase variation. When we write 7, in the 
approximate form r,=2,+ PS eqn. (2) becomes 


Been ty =r [ {| Fe yexp 4) + ym) 
x Ge ee ies iat). ear ee a 
ma, r 
yee oe , (8) become, respectively, 
BR ee NF ETE ari exp (Sa (+n?) dé dr, 
(15) 
Bet." _ FF (16) 
Cree Say eel eee 
where 
pe, FA 
eis, 
In the integrals A and Bb, zis replaced by z’. The eqns. (11) become 
fee QO = E.= a F-. 
pe eee ie 
& Sere a iT 3 Ty 2 a) 2 
ee Ge oe 1) A ei a eas) | 


3. THE JOINT GAUSSIAN PROBABILITY DISTRIBUTION 


One may show by calculating the nth moment of #, and £, that although 
the variances of H, and H, change considerably ure the course of 
propagation, the Fores of ine probability distribution functions changes 
only a little. This result will be presented in detail elsewhere. 

If F and /’, have a joint Gaussian distribution we are therefore approxi- 
mately: justified in assuming that #, and H, also have a joint Gaussian 
distribution. 
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ee Se 2 
f(y, Bs) eo “wo 


eee tea aCe AUP tear 19) 
2 


— 
e 
4 
pb 
Dp 


where 


O 


OA is the average signal EB. The line BC is the ‘ ellipse ’ f=const. in the case 
of amplitude modulation ; the circle DEF is the ellipse for a randomly 
phased signal. 


f =const. is an ellipse in the #,, #, plane, with centre at (Z,, #,). Ampli- 
tude modulation is indicated by a straight line parallel to the average 
signal. The fading signal is said to be randomly phased when / = const. is 
a circle. See fig. 1. 

If we are to discuss effectively the way in which # becomes randomly 
phased as the wave propagates, we require a measure of the ellipticity 
which is independent of the particular choice of axes H,, H,. We cannot 


P.M. 3D 
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use + because when the curve f= const. is not circular r can always be made 
to vanish by a suitable rotation of the axes. 7=0 is a necessary but not 
sufficient condition that f =const. be a circle. We therefore define an 
intrinsic correlation py of the fading signal as follows: py is the maximum 
value of 7 under rotation of the axes H,, E,. 

The problem is to find py in terms of o,, o2, 7. We will first transform the 
coordinates H,, H, to a system A,, A, whose origin is at the centre of the 
ellipse and whose axes lie along the major and minor axes of the ellipse. 
Let the transformation be as follows. See fig. 2. 


E,-E£,=), cos —A, sin 6, } (19) 


HE, — EH, =), sin 6+A, cos 6. 


Fig. 2 


O 


OA is the average signal H. The axes )j, 2 are along the axes of the ellipse. 
The correlation coefficient is maximum in the py, uw, coordinate system. 


When substitution is made in the quadratic form in (18) we obtain 


12 (= 7] in sin? @ 27 sin 6 cos *) Ae (= 6 re cos? @ e 27 sin @ cos *) 


1 2 2 
O71 O»5 0109 


a 


2 2 
O71 Oy 0109 


= 2 2 


é cos@sin@ cos @sin @ T 
re DAS Ng a 
O71 om 010 


(cos? 6 — sin? a). (20) 
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The axes A,A, lie along the axes of the ellipse when the coefficient of ),A, 
in (20) vanishes. Therefore 


270,0 
tan 20=—-— 42 ww. QD) 


go 
The frequency function of ,, A, is 


Fy de) = <—Kesoxp { 4 ‘e soll. al be) 
Qar/ Ayr," 1 ro? 


A further rotation will be made to determine the system of coordinates in 
which the correlation coefficient 7 isa maximum. See fig. 2 


A, =}, COS $— py sin 4, \ 
A=, Sin J + py Cos d. 

The correlation coefficient of 14, 4, 7, say, is obtained from the coefficient 

of p42 In the quadratic form 

Ay? Ag? > (cos? sin? (ee cos? 

AZ OAS Whe 2 do? 


—2 a = us sin d COs p 
Mike AZ 3 re? 3 


(23) 


whence 

a cos ¢ sin o(1/Ay?— 1/A,?) 

~ {1/(AP A) + sin® f cos? $(1/A.2 — 1/2g2)222 
7, 18 Maximum when ¢=7/4, clearly. one 


2 See? 
= AAPA? ef Tee es (24) 
LEAF Ay A +e 
This shows that if « is the eccentricity of the ellipse, then |p,| = ¢?/(2—?). 
The ratio ),2/A,? is the ratio of the coefficients of A,? and A,? in (20). After 


some simplification we obtain 


a 


Py = Max T 


Dye pg 1 
— O71 29— {( Pe 22+ (2 0.0 yep, (25) 
= —— sec 240 = ——_ 71 (9. O95 iP 192 
PN Ge +o? ee ag 2 


As a last point we give the expression for /(A,, A 


2) 
1 [ 1 (Z Ay ree nee —)\ (26) 
A m1(o4” + o57)(1 — py”) ee oy +o,*\l+py 1—py 
The propagation of the fading amplitude will be described in terms of 


py and 6. 


in terms of o,, 09, py, 


4. Tue InTRINSIC CORRELATION OF THE FADING AMPLITUDE 
py and @ as defined by eqns. (21), (25) are functions of z, since H, and £, 
neperd on z. We shall find p,y(z) and 4(z) in terms of the variances of 
F,F, and the A’s and B’s. Let py(0), 6(0) denote the corresponding 
Peepectics of the ellipse of F,F,. We are free to choose the orientation 


2D2 
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of the ellipse of /',F', and shall suppose that the axes FF, are parallel to 
those of the ellipse. Then from eqns. (11), (18), (21), (25), 


2 2 rs 
tan 20(z)= — ee =- = (Hii, —* 2) 
Caer (HP Het Ee?) 
_ —(FP-F)B, + (FP -F,?) Bs 27) 
— (Fy? —F,)A, + (f° —F,7)A, 
Be. N\ (hee 
py(2) os 71 sec 20= oe = = =i sec 26 
Ce pase (Hp? —H?) + (Hy? — #,") 
ies eae at sec 26(z). (28) 


For a point source it is only necessary to nite zin the A’s and B’s by z’. 
These results are much simpler when the three spatial correlation 
functions of F, and F’, are identical. This is the case when the same 
diffraction mechanism, such as a thin screen of varying refractive index, is 
responsible for the fluctuations in both fF, and F,. If two diffracting 
layers modified the amplitude and phase separately the spatial correlations 
would be different. Suppose then that A,=A,=A, B,=B,=B. We 

then have 
tan 20= — B/A, 5) a een 


(Fy? —F,?) —(F2—F,) 
(FF?) + (F2—F,”) 
= py(0)/(A2 + B?). Ee ne ts bee (Ei), 


py(0) is defined to be the intrinsic correlation of the fluctuating amplitude 
F’, which is seen from (24) to be 


jeg A sec 20 


OP Pee ey) 
eae ieee) 
The +root is chosen in (30) so that py(z)>py(0) as z+0. Equations 
(29), (30) are independent of the system of axes used for # and F' provided 
@ is the angle between the major axes of the ellipses. 
These equations have a simple physical interpretation. According to 


(10), 


—A+i1B= ml}. pl, Jen (SE 


i ec plv/2E, v2nhexp (SZ (e+ x) dé dr 
(31) 


Therefore —A+iB is the complex amplitude at the centre of the Fresnel 
diffraction pattern by an irregularity ‘over which the amplitude is 
p(/2&, 2m) . py(2)/py(0) and 6(z) are therefore the amplitude and phase 


t (g2 +7) dé dn 
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of the diffraction pattern. The variation with distance z of the correlation 
at the receiver is a sort of focusing in the diffraction. The fading is randomly 
phased when py(z)=0, i.e. when the signal at the receiver due to a single 
irregularity has vanished. 

Let us consider the previous example (eqn. (12)) 


1 
pi=pe=pa=exp| ~ SG +9%) |. 
; 0 
Equations (29), (30) become 


L a ; 4)222\ —1/2 
Py(2)= Vata x)= (: ot ra) pPy(), | 
Lay ee ft 
we To” | 


(32) 


As 2 00 the ellipse of H becomes circular and rotates through an angle 
of 7/4. 
Fig. 3 


A 


O 


This shows the construction of the modulation imposed on the signal at the 
maximum height. The modulation is the line AB, the diagonal of the 
ellipse. The ellipse rotates counter-clockwise as z increases. 


From the definition of py(z) we know that |py|<1. This implies a certain 
maximum value for 2, viz. 


TT 2 1 
max — : 2 Fue 
2A py’ (2) 
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We therefore see that if the irregularities are known in size and there is a 
known correlation in the fading at the receiver, then there is a maximum 
height at which the irregularities can be supposed to be situated. At that 
maximum height py(0)=1. 

There is a simple geometrical relation between the straight modulation 
at the maximum and the ellipse at any other distance. From (32) we 
deduce that p,y(z)=py(0) cos 20. Therefore the straight line is the diagonal 
of the ellipse as shown in fig. 3. 


4.1. Effect of Point Source 


When a wave from a point source is incident on the screen then (27) and 
(28), etc. are still valid provided one changes z to z’ (see (16)). 


Fig. 4 


Typical daytime fading record on 16 ke/s : 10.00-11.00 G.M.T. July 6, 1949. 
The points represent the signal at minute intervals over the period of an 
hour. The diagonal of the probability ellipse is drawn to show the phase 
modulation at the maximum height. 


4.2. Hxperimental Results 


A certain amount of experimental work has been done on 16 ke/s. 
The signal from GBR, at Rugby, is recorded after nearly vertical reflection. 
The receiver records the components of the sky-wave which are in-phase 
and in-quadrature with the ground wave. By plotting these components 
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as Cartesian coordinates one is able to study the fading signal in the light 
of the foregoing analysis. In fig. 4 an example is shown. The points 
represent the signal at minute intervals over the period of an hour. The 
curves of constant probability density are nearly elliptical. When the 
diagonal of the ellipse is drawn it appears that the modulation at the 
hypothetical maximum height is, in this example, very nearly phase 
modulation. 

It is difficult to make a survey of the numerous records hitherto obtained 
because most of them show a large systematic phase variation which 
obscures the more rapid short period fading. The fading has a period of 
about three minutes, while the slow phase change may amount to 7/4 in 
an hour. 

The whole theory is at the limit of validity in the case of 16 ke/s, since 
r=5 A, and eqn. (1) is true under the assumption that r>A. However the 
only extensive data are for this frequency and it is expected that the result 
will be not far from the truth. 

Five hourly periods have been analysed. Four of these, including the 
one in fig. 4, show very nearly phase modulation. The other appears to 
be a mixture of phase and amplitude modulation. 

In the example shown the parameter x may be calculated from the 
angle 6. From the relation py(z)=py(0) cos 20 we have 20=46°. There- 
fore, from (32) «=cot 26=0-96. Assuming 2’ =z/2=45kmandA=18-8km., 
we obtain r9=23 km. This is in tolerable agreement with the estimates 
made by the spaced receiver method of measuring irregularity size 
(Hargreaves, private communication). 


5. CONCLUSION 


The fluctuating in-phase and quadrature components of a fading signal 
have been assumed to have a joint Gaussian probability distribution. 
This distribution has not been described in terms of the usual variances and 
correlation coefficient, but in terms of two parameters py and 6: py 
measures the extent to which the signal is randomly phased and 6 measures 
the change in phase of the modulation on the signal after it has left the 
screen. These two parameters are the amplitude and phase of the Fresnel 
diffraction pattern due to a single irregularity. The experimental records 
of 16 ke/s fading analysed according to this theory show that at this 
frequency the modulation is highly correlated at the ionosphere and that 
it is generally phase modulation. 
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CORRESPONDENCE 


The Hardening of Copper by Neutron Irradiation 


By M. A. Apams and P. R. B. Higerns 


Metallurgy Division, Atomic Energy Research Establishment, Harwelt 
[Received February 9, 1959] 


Varrous studies have been made of the effects of neutron irradiation 
on the mechanical properties of metallic materials (see the review articles 
of Cottrell (1956, 1957) Glen (1955) and Kinchin and Pease (1955)), and 
it is now generally recognized that irradiation may cause hardening in two 
broadly different ways: (a) by making the nucleation of slip difficult 
(source hardening) and/or (b) by making slip propagation difficult (lattice 
friction hardening). Copper has been most extensively investigated but 
in none of the experiments on this material has the relative contribution 
of the two forms of hardening been clearly separated. One means of 
making the separation is, however, suggested by work that has been done 
on mild steel, for which it has been shown, both theoretically (Petch 1953) 
(Cottrell 1958) and experimentally (Hall 1951) (Petch 1953) (Hull and 
Mogford 1958), that the lower yield stress is related to the grain size by 
the equation : 

o,=o;,t+k,.d-'? Pieper er ag ees MES 


where, following Cottrell’s (1958) notation, d is the half-grain diameter, 
co, the shear stress at the lower yield point, k, a measure of the shear 
stress to operate dislocation sources and o; the shear stress resisting the 
movement of dislocations across the slip plane. Equation (1) describes 
generally the conditions for yield propagation when there is a sharp 
yield point. Copper shows such a yield point after irradiation (fig. 1.) 
so that measurements of the yield strength as a function of d~"? on speci- 
mens irradiated various amounts should reveal whether source hardening 
(changes in /,,) or lattice friction hardening (changes in o;) occurs. 

Wires of 0:085in. diameter were prepared from Matthey spectrographic 
(99-999°%) copper rods by successively swaging (50% reduction in area) 
and vacuum annealing (600°c for 2 hours). After a final swage to 0-060 in. 
diameter, batches of 3-5cm long specimens were vacuum annealed for 
2 hours at temperatures from 400°c to 1000°c to give grain sizes (deter- 
mined from a metallographic examination of the specimens after testing) 
in the range from 0:0025cm diameter to 0-08cm diameter. Specimens 
covering the complete grain-size range were sealed in evacuated silica 
capsules and irradiated at a temperature of about 50°c for 1 and 3 weeks 
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Lower yield stress values at 20°c and —78°c of copper samples covering a 
range of grain sizes, after irradiations of 1 and 3 weeks in a fast neutron 
flux of 3 x 1011 n cm- sec-!. 
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inside a fuel element in the sero reactor. The measured thermal 
neutron flux was 4:5 x 104 nem-?sec—!, the measured flux above about 
Imev, 3x 10"nem~ sec", and the estimated flux between 0-4 ev and 
1 Mev, 1-5 x 10!2n cm? sec (Wright 1959, private communication). 

Tensile tests were made at 20°c and —78°c in the autographically 
recording tensile machine described by Adams (1958). The wires were 
held by their ends in carefully aligned clamp type grips leaving a gauge 
length of 2cm, and were extended at a strain-rate of 5-7 x 10-4see~!, 

Results for the two irradiations and two testing temperatures are 
shown in fig. 2 as a plot of lower yield stress against d-¥? and can be 
summarized as follows: 


1. At a fixed neutron dose and testing temperature the yield stress of 
irradiated polycrystalline copper increases, within the limits of 
accuracy of the experimental measurements, linearly with increasing 
d-? in agreement with egn. (1). 


2. The effect of increasing irradiation is to raise both o, and k,, (compare 
curves | and 3 and curves 2 and 4). Thus lattice friction hardening 
and source hardening proceed simultaneously. 


3. The increased yield strength of irradiated copper at — 78°c compared 
with 20°c, arises from a temperature dependence of the lattice friction 
component of hardening (compare curves 1 and 2 and curves 3 
and 4). 


The feasibility of the various theories of radiation hardening can now 
be assessed in the light of these results. Thus, it seems unlikely that a 
source hardening mechanism with the required lack of sensitivity to 
temperature can be provided by a point defect atmosphere locking effect, 
or by a simple jog process. If, however, dislocations become heavily 
jogged, as recent observations on quenched metals suggest (Hirsch et al. 
1958), temperature independent locking is conceivable, since a large length 
of jog on some plane other than a close-packed one may constitute an 
almost immovable lock. Vacancies condensed as voids on dislocations 
(Coulomb and Friedel 1957) could also cause source-hardening of a type 
which is relatively temperature insensitive. Regarding the lattice 
friction component of hardening, the ‘peg’ theory (Cottrell 1957) seems 
most likely to provide a mechanism with the required sensitivity to 
temperature. Aggregates of defects (Kunz and Holden 1954a, b). 
dislocation rings formed by the collapse of vacancy clusters (Smallman 
and Westmacott 1959), and disturbed zones in the lattice (Seeger 1958) 
are possible ‘pegs’ through which gliding dislocations could be forced 
by the applied stress with appreciable assistance from thermal fluctuations 
at relatively low temperatures. 

Further work is proceeding to determine the dose dependence and 
annealing characteristics of the two types of hardening. Full details of 
the experiments will be published later. 
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Variation of Young’s Modulus of Fused Quartz Fibre with Diameter 


By D. J. Trirron 
Cavendish Laboratory, Cambridge 


[Received April 10, 1959] 


THERE has been some dispute (Jurkov 1935, Reinkober 1939, Tighe 1956) 
as to whether Young’s modulus of fused quartz-fibre varies with diameter. 
Some recent experiments by the author demonstrate clearly that it does, 
though they do not fill in the detailed form of the variation. 

The observations were made during the development of a quartz-fibre 
anemometer for measuring the low air speeds of free convective flows. 
This instrument consists of a fused quartz-fibre (of diameter in the range 
20 to 100 microns and length in the range 1 to 3cm) cemented into a 
piece of hypodermic tubing. When this is placed in an air flow the 
deflection of the free end, observed on the eyepiece scale of a tele-micro- 
scope, is a measure of the velocity. 

During calibration runs, it became clear that the results could not be 
consistently interpreted using the value of 5-18 x 1044dynesem- for 
Young’s modulus of quartz-fibre in Kaye and Laby’s (1956) tables. Hence, 
some measurements were carried out. The method was a determination 
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of the resonant frequencies of some of the fibres, as mounted for velocity 
measurements. A fibre may be regarded as a bending beam with one 
end free and the other clamped. Hence, its equation of motion is 


oY ory 


which becomes 


4 
EI oo —47*mvy, = Fy 

for a periodic motion of frequency v. Here y is the displacement of the 

fibre at distance « from its fixed end and y, the amplitude of this displace- 

ment; F' is the bending force per unit length and F, its amplitude; m, 

and J are the mass per unit length, Young’s modulus, and moment of the 

cross section of the fibre. 'The boundary conditions on this are 


dy 
{} ie Se =—() i =i). 
Yo . at & 
Byy _ BY o 
eee (ety =| 
da? das nae ee 


so that the resonant frequencies are given by 


il 1] i i 

Waa : ’ . . =() d= (eee | 
e? e? —e? —¢@ ? EI 

e? =—¢* —iet ie 


which reduces to 
cos ¢cosh¢= — 1. 
The first four eigenvalues of ¢ are 
$,=1-875, 6,=4:712, $3 =7-853, $4= 11-00. 
Young’s modulus is thus given in terms of the resonant frequencies as 
"f 4777l4mv,,? a 6477I4p9 [vy \? 
Tone a? One 

where py is the density of quartz and d the diameter of the fibre; v, is the 
observed nth resonant frequency. This treatment ignores damping, but 
the fact that good agreement was obtained in the experimental observa- 
tions of different eigenvalues is a justification for neglecting its effect on 
the resonant frequencies. | 

The resonances were found by mounting the fibre 10cm or so in front 
of a loudspeaker connected to an electronic oscillator. The fibre some- 
times responded to one of the higher harmonics in the loudspeaker’s note ; 
it was then necessary to discover at what frequency the fibre was vibrating, 
which was done by observing it in the light of a stroboscope. The fre- 
quencies were all in the low audio range, and there was overlap between 
the ranges covered by the first few eigenvalues for the various fibres ; 
the results are thus not due to some frequency-dependence. 
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Of the other quantities needed for calculating H, the diameters of the 
larger fibres were measured with a travelling microscope and those of the 
small ones by a wedge-fringe optical interference method. 

The fibres investigated came from two different sources. Ones of 
approximately 20 and 45 microns in diameter were supplied by Messrs. 
Taylor, Taylor and Hobson Ltd. of Leeds who kindly measured the density 
of their original material for me (2-20gem~*). Since the experiments of 
Douglas and Isard (1951) showed that density variations produced by 
heat treatment are only of the order of 0:2%, this figure may be used for 
the fibres. A fibre of approximately 75 microns in diameter was made from 
some scrap fused quartz, of which the origin is not known. ‘The same 
figure was used for the density of this; Tighe (1956) remarks that there 
are no significant variations in the density of different samples of fused 
quartz. 

The results were 


Fibre diameter Young’s modulus 
(cm) (dynes cm~?) 
por x 108 8-1 104 
4-42 x 10-3 6-85 x 101 
oe x 10 5-95 x 101 


I estimate that these values for Young’s modulus are accurate to around 
4%. Each is based on the results for a single fibre, but, further experi- 
ments were made with fibres of other lengths from the same reels. 
Accurate diameter measurements were not made on these (since this is a 
tedious and time-consuming job), but, supposing them to have the same 
diameters gives agreement with the above results to within 3%. It is 
thus safe to conclude that the above are not freak results produced by, for 
example, a flaw in one of the fibres. 

It is doubtful whether samples of quartz from different sources always 
have the same constitution ; Drane (1929), for instance, concluded that the 
results of various experimenters on the elasticity of fused quartz rod could 
be explained only by actual variability. Hence, a comparison between the 
75-3 micron fibre and the others is of doubtful value. On the other hand, 
the manufacturers assure me that the 19-7 and 44:2 micron fibres were 
made of the same material, so that the difference, considerably more than 
the experimental error, in Young’s moduli of the two is the most interesting 
result of this work. The difference is very likely due to the rate of cooling 
of the fibre during manufacture. Comparisons with other workers are 
also open to the criticism that the material may not be just the same. 
It is, however, perhaps worth noting that my values of H for the 44:2 
and 75-3 micron fibres are within the range of the various curves given 
by Reinkober (1932, 1937, 1939) but that the increase as the diameter is 
decreased is even larger than that observed by him. Jurkov (1935) 
observed a modulus of about 9-4 x 10!dynesem- for all diameters; 
this seems unlikely in view of my results. 
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The resonance method can be recommended as a means of finding the 
full form of the variation. However, it should be noted that if there is a 
variation in the local elasticity of an individual fibre with distance from the 
axis the modulus would be different for bending than for stretching. 

The results of this work have been successfully used in determining from 
calibrations of a number of quartz-fibre anemometers the drag produced 
by low Reynolds number flow past cylinders. A report on this will be 
published elsewhere shortly (Tritton 1959). 
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' Waves in Electron Clouds 


By J. R. Pierce and L. R. WALKER 


Bell Telephone Laboratories Inc., Murray Hill, 
New Jersey, U.S.A. 


[Received February 2, 1959] 


In a recent paper in this journal (Piddington 1958) the author in referring 
to an earlier paper of ours (Pierce and Walker 1956) quotes us as saying: 
‘ No definitive conclusion can be drawn from it (the steady state analysis) 
alone.” What we said was: ‘“‘ This relation, the dispersion relation, is 
obviously significant for the problem on hand since it describes an intrinsic 
property of the ion stream and may, of itself, give some valuable infor- 
mation. But no definitive conclusion can be drawn from it alone about 
any specific problem, since only the boundary conditions can say what 
part of the wT spectrum is excited and to what extent.’ 
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REVIEWS OF BOOKS 


Modern Geometrical Optics. By Max Herzpercer. (New York and London: 
Interscience Publishers, 1958.) [Pp. 504.) $15. 


Tuts book, by one of America’s leading optical specialists, deserves an honoured 
place in every optical library. It contains thirty-nine chapters (four in an 
Appendix), a bibliography and an index. In his preface, the author reminds 
us that lens design is still more of an art than a science and that the acquaintance 
of the average lens designer with higher mathematics is limited. For this 
reason, he has tried to restrict the demands on the reader’s mathematical 
technique to the minimum. But this minimum is not a very small one and the 
average lens designer who reads the whole of Dr. Herzberger’s book will find 
that he has learned a good deal of higher mathematics as well as of geometrical 
optics. Without this mathematics, a systematic account of the subject would 
be impracticable. Vectors are freely used, for the sake of the simplicity which 
they bring, and an admirably lucid introductory course on vector analysis is 
provided in the Appendix by Dr. Erich Marchand. Another chapter in this 
Appendix, entitled Miscellaneous Mathematical Tools, supplies necessary back- 
ground knowledge of matrices, of the method of least squares, of Gaussian 
brackets, and of polynomial approximation methods. By thus allowing the 
reader to educate himself as he goes along, Dr. Herzberger is able to proceed 
from the thoroughly practical and mathematically elementary ray-tracing 
methods which occupy his first seven chapters, through Gaussian optics in 
the next seven, to a concise and elegant presentation of Hamilton’s ideas, 
according to which the laws of geometrical image formation by an optical 
system are expressed through a single characteristic function, or Eikonal 
function, as it later came to be called. 

Of many other interesting topics treated in the last part of the book we here 
mention the method of fitting a fifth-order mathematical model (simplified 
approximation) to an optical system on the basis of the results of tracing a 
relatively small number of meridonal and skew rays, and the use of spot 
diagrams as a tool in lens design. The author’s use of diapoints often lends an 
agreeably modern flavour even to the more elementary discussions. 

Dr. Herzberger’s wide knowledge and enthusiasm for his subject illuminates 
the whole book and in the fourteen years spent, as he tells us, on its writing 
and rewriting he has produced a major contribution to the optical literature. 

EK. H. Lryroot. 


Theories of Figures of Celestial Bodies. By WENCESLAS 8. JARDETZKY. (New 
York: Interscience Publishers, 1959.) [Pp. 186.] 45s. 
A MATHEMATICIAN who wants to find out what has been done, and who has been 
doing it, in the subject of the mechanics of rotating fluid masses, could hardly 
do better than start by consulting this book. . 
Professor Jardetzky’s breadth of knowledge of the European and American 
literature on the subject is evident in this concise and clearly written historical 
review. In the first part of the book, the author gives an outline of the work 
done by early investigators, Jacobi, Poincare, and many others, on bodies of 
fluid rotating as if solid. The basic mathematical equations are stated, methods 
and results of solutions are indicated, physical principles and various simplifying 
assumptions are discussed. In the second part, more complicated problems are 
considered, particularly such conditions as are known to exist in actual celestial 
bodies—zonal differences of rotation and other internal motions, bodies of mixed 
composition, small oscillations, and the effects of external fields of force. 
There is a very extensive bibliography. J.E. J. 
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D. HULL and D. E. RIMMER Phil. Mag. Ser. 8, Vol. 4, Pl. 77. 


Voids along grain boundaries in a specimen fractured at 480°C with a uniaxial 
stress of 4730 p.s.i. Hand polished and etched. x 270. 


sins 


x 2000. 


Same as fig. 3. 


BD, HULL andeDaEs RIMMER Phil. Mag. Ser. 8, Vol. 4, Pl. 78. 


Fig. 5 


Voids formed at 410°c with stress of 4730 p.s.i. x 270. 


Fig. 6 


Voids formed at 410°c with stress of 3000 [Ossie 3 Bri). 


A. E. BROWN and E. G. RICHARDSON Phil. Mag. Ser. 8, Vol. 4, Pl. 79. 


Fig. 8 


(c) (d) 


Schlieren pictures of mixture at various temperatures : 


(a) 3129 “(b) 32°; (c) 33°: (d) 36°e. 


P. B. HIRSCH, et al. Phil. Mag. Ser. 8, Vol. 4, PI. 80. 


0-5 u 
| 
Contrast produced by surface slip steps. At A and B dislocations are 
associated with slip steps. Mag. x 50 000. 


P. B. HIRSCH, et al. Phil. Mag. Ser. 8, Vol. 4, Pl. 81. 


0-5 
—— 
Dislocation networks at the surface. At A dislocations meet the protected 
surface in visible lines corresponding to slip steps. Mag. x 50 000. 


Surface of a fatigue specimen. Note at A offsets produced by slip on two 


systems, at B slip similar to that in simple tension, at C fatigue striation 
and at D patches of stacking fault. Mag. x 33 000. 


Py Bw hIRSCH 6 al: Phil. Mag. Ser. 8, Vol. 4, Pl. 82. 


Slip steps at the edge produced by dislocations lying in the slip planes. Note 
shifts in extinction contour at A and B. Mag. x 45 000. 


Pao 
HIRSCH, et al. Phil. Mag. Ser. 8, Vol. 4, Pl. 83. 
Fig. 4 


Intensity change on crossing slip steps. Note also the offsets produced where 
the slip lines intersect. These can be compared with the steps at the 


edge. Mag. x 55 000. 


P. By SIRSC;. et ai. Phil. Mag. Ser. 8, Vol. 4, Pl. 84. 


Hig. 7 


Fatigue specimen showing at A slip line without associated edge steps and at 
B light regions corresponding possibly to fatigue cracks. Mag. x 45 000. 


Gao riIRSCH, et’ al, Phil. Mag. Ser. 8, Vol. 4, Pl. 85. 


Fig. 8 


0-5 
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i iatl Al osely spac shi Arrow shows direction of 
Fatigue striation showing closely spaced slip. ae 
: maximum taper magnification. Mag. x 45 000. 


E20. MALE Phil. Mag. Ser. 8, Vol. 4, Pl. 86. 
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»MnAl, annealed at 400°c for 7 days ( x 550). 


E/Os HALL Phil. Mag. Ser. 8, Vol. 4, Pl. 87. 


Ag,MnAl, annealed at 400°c for 7 days ( x 550). 


Fig. 8 
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Ag,MnAl, annealed at 400° for 7 days ( x 550). 


E. O. HALL Phil. Mag. Ser. 8, Vol. 4, Pl. 88. 


Ag,MnAl, annealed at 400°c for 7 days ( x 550). 


E. O. HALL Phil. Mag. Ser. 8, Vol. 4, Pl. 89. 


Ag;MnAl, quenched and aged at 310°c for 7 days ( x 2000). 
(Slightly defocused.) 


Fig. 12 


Ag;MnAl, quenched and aged at 400°c for 7 days (x 2000). 


E; O) HALL Phil. Mag. Ser. 8, Vol. 4, Pl. 90. 


Fig. 13 


Ag;MnAl, quenched and aged at 240°c for 7 days. Carbon replica ( x 10 000). 


Fig. 14 


Klectron diffraction pattern from same specimen. 


E, O. HALL Phil. Mag. Ser. 8, Vol. 4, PI. 91. 


Fig. 15 


Ag;MnAl, quenched and aged at 400°c for 7 days. Carbon replica 
: (x 10 000). 


